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Abstract. The purpose of this work is to resume investigation of Galactic thin disk dating using nucleocosmochronology with 

Th/Eu stellar abundance ratios, a theme absent from the literature since 1990. A stellar sample of 20 disk dwarfs/subgiants of F5 

to G8 spectral types with -0.8 < [Fe/H] < +0.3 was selected. In stars with such spectral types and luminosity classes, spectral 

synthesis techniques must be employed if we wish to achieve acceptably accurate results. An homogeneous, self-consistent 

] set of atmospheric parameters was determined. Effective temperatures were determined from photometric calibrations and 

Ha profile fitting; surface gravities were obtained from Teff, stellar masses and luminosities; microturbulence velocities and 

metallicities were obtained from detailed, differential spectroscopic analysis, relative to the Sun, using equivalent widths of Fei 

, and Fell lines. Chemical abundances of the elements that contaminate the Th and Eu spectral regions (Ti, V, Cr, Mn, Co, Ni, 

' Ce, Nd, and Sm) were determined through spectroscopic analysis. Abundance uncertainties were thoroughly scrutinised, their 

average value - (0.10 ± 0.02) dex - being found to be satisfactorily low. Eu and Th abundances were determined by spectral 

^ ^, synthesis of one Euii line (4129.72 A) and one Thii line (4019.13 A), taking into account the detailed hyperfine structures 

I ' of contaminating Co lines, as well as the hyperfine structure and isotope shift of the Eu line. Comparison of our abundances 

O I with literature data shows that our results exhibit a similar behaviour, but a considerably lower scatter (36% lower for Eu, and 

61% lower for Th). The [Th/Eu] abundance ratios thus obtained were used, in the second paper of this series, to determine the 

^ ' age of the Galactic disk. 

■ 
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1. Introduction jects and events. The Th/Eu chronometer was first proposed by 



je, 



( 1989) as a way of assessing the age of the Galactic disk. 

Modern methods of estimatmg the age of the Galactic thin ,^ • j ^ . 2i2rru u t a nc u iri-r 

, o o main advantages are that Th has a 14.05 Gyr halr-Me, 

disk , like dating the oldest open clusters by isochrone fitting . c j c j f.u u • j j 

° ^ -' ° I.e., or the order or magnitude or the age being assessed, and 

and whites dwarfs by cooling sequences, rely heavily on stel- . t- j f . i • u ■ 

-' to 1 'J J lYiat Eu provides a satisfactory element for compa rison, being 

lar evolution calculations. Nucleocosmochronology is a dat- j j i , i • i /nnm j- . 1^ ■ — TTl 

produ ced almost exclusively (97%, according to IBurris et alJ 

ing method that makes use of only a few results of main se- |L„„4, , , . j n 

■' I2()()(ll) by the same nucleosynthetic process that produces all 

quence stellar evolution models, therefore allowing a quasi- rn. j . . / x • j- 

^ o 1 Yjj^ jjjg j-apj(j neutron-capture process (r-process). Its main dis- 

independent verification of the afore-mentioned techniques. j , i - • j vcc » i. . • . 

^ ^ advantage lies in the difficulty one encounters when trying to 

Nucleocosmochronology employs the abundances of ra- determine stellai- Th abundances. This difficulty arises from the 

dioactive nuclides to deter mine timescales for asti-ophysical ob- equivalent width (EW) of the Th ii line used, and from 

Send offprint requests to: E.F. del Peloso the fact that it is severely blended with other much stronger 

* Based on observations collected at the European Southern lines. 
Observatory, La Silla, Chile, under the ESO programs and the ESO- 

Observatorio Nacional, Brazil, agreement, and at the Observatorio The first to suggest the presence ofaThiiUne in thesolar 

do Pico dos Dias, operated by the Laboratorio Nacional de spectruin located at 4019.137 A, werelSitterly & King|(|l943|). 

Astroflsica/MCT, Brazil. ISevernvl ( Il958l) was the first to measure the EW of this line, de- 

' All references to the Galactic disk must be regarded, in this work, riving an upper limit for the solar thorium abundance. The first 

as references to the thin disk, unless otherwise specified. confirmed Th detection and abundance determination in a star 
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other than the Sun was accomplished bv lCowlev et alJ d 19751) 
in HR 465, an Ap star. Until today, only two works available 
in the literature have pre sented Th abundances for samples 
of Galactic disk stars: da Silva et al.' ( 1990) and iMorell et alJ 
(1992, MKB92). Of these, only da Silva et al. present a chrono- 
logical analysis. 

This project aims at resuming investigation of Galactic disk 
dating using [Th/Eu] abundance ratios, a theme absent from 
the literature since 1990.^ In this paper. Part 1 of a series, we 
carried out the preliminary, observationally-oriented steps re- 
quired. First, an appropriate stellar sample had to be elected. 
Selection considered the suitability of the sample to the task of 
Galactic disk dating, with criteria like metallicity and spectral 
type range. In stars with spectral types and luminosity classes 
adequate to our study, which are F5-K2 dwarfs and subgiants 
(see Sect. 12. Il l, there are only one Th and one Eu line adequate 
for abundance determinations. This means that spectral syn- 
thesis techniques had to be called upon so as to obtain accept- 
ably accurate results. As prerequisites to the synthesis, we de- 
termined accurate atmospheric parameters and chemical abun- 
dances of the elements that contaminate the Th and Eu spectral 
regions (Ti, V, Cr, Mn, Co, Ni, Ce, Nd, and Sm). Eu and Th 
abundances were determined for all sample stars by spectral 
synthesis, using our atmospheric parameters and abundances. 
The [Th/Eu] abundance ratios thus o btained were used, in the 
second paper of this series ( del Pelos o et alJEo05a. Paper II), 
to determine the age of the Galactic disk. 

2. Sample selection, observations and data 
reduction 

The st ellar s ample was define d using the Hipparcos cata- 
logue jPerrvm an &F,SaI |1997I) . upon which a series of se- 
lection criteria were applied. All objects in the final sam- 
ple were observed with the Fiber-fed Exten ded Range Optical 
Spectrograph (FEROS; iKauferet al.ll 19991) fed by the 1.52 m 
European Southern Observatory (ESO) telescope, in the ESO- 
Observatorio Nacional, Brazil, agreement. Spectra were also 
obtained with a coude spectrograph fed by the 1.60 m tele- 
scope of the Observatorio do Pico dos Dias (OPD), LNA/MCT, 
Brazil, and with the Coude Echelle Spectrometer (CES) fiber- 
fed by ESO's 3.60 m telescope and Coude Auxiliary Telescope 
(CAT). The selection criteria are related in detail below, fol- 
lowed by the description of the observations and their reduc- 
tion. 

2.1. Selection criteria 

Multiple criteria were applied to the Hipparcos catalogue, com- 
posed of 118 218 objects. Initially, we eliminated objects with 
parallaxes lower than 0.010", in order to ensure the minimum 
uncertainty in the derived bolometric magnitudes (22 982 ob- 

^ In this paper we obey the following customary spectroscopic no- 
tations: absolute abundance log£(A) = log,Q(A'A/A'H) + 12.0, and 
abundance ratio [A/B] = logi(,(A'A/A'B)star - logioCA'A/A's)©, where 
A'a and are the abundances of elements A and B, respectively, in 
atoms cm"^. 



jects left). Only stars with declination <5 < -i- 20° were kept, to 
allow the observations to be caiTied out in the southern hemi- 
sphere (15 898 objects left). Stars fainter than visual magnitude 
V - 7.0 were removed, so that high resolution, high S/N ratio 
spectra could be obtained with relatively short total exposures 
(i.e., less than 3 hours in total) on small and medium sized tele- 
scopes (3272 objects left). 

The subsequent criteria aimed at constructing the most 
suitable sample for Galactic disk nucleocosmochronology. 
Spectral types were restricted to the F5-K2 range because 
these stars have life-times comparable to the age of the Galaxy 
(1744 objects left). Only luminosity classes IV and V were al- 
lowed, to avoid stars whose photospheric abundances had been 
altered by dredge-up episodes, and to minimise non-LTE ef- 
fects (925 objects left). Since spectral analysis would be dif- 
ferential, we limited the sample to stars with atmospheric pa- 
rameters similar to the Sun by selecting objects with a (B-V) 
color index in the interval [+0.45, +0.82], so that their effec- 
tive temperatures fell in the range Tetfo + 400 (744 objects left). 
As the ultimate objective of this work is the determination of 
the age of the Galactic disk, we eliminated stars with [Fe/H]< 
-1.00 (252 objects left); for this purpose, we employed a ver- 
age rn etallicities from the catalogue of Cayrel de Strobel et alJ 
il200ll) . 

Stars listed as double in the Bright Star dWarren & HoffleiJ 
ll287) or Hipparcos catalogues were rejected. As a last cri- 
terium, we only kept stars whose masses could be determined 
with the Geneva sets of evolutionary tracks ( Schaller e t alJ 
1292, Charbonn el et al 1993, Schaerer et al. 1993a b, and 
Charbonnel et alJ 1 19961 hereafter collectively referred to as 
Gen92/96). For this purpose, stars were required to be located 
between tracks in at least two of the HR diagrams constructed 
for different metallicities. This left 157 objects in the sample. 

To reach a suitably sized sample, we kept only the brightest 
stars to end up with 2 to 4 objects per 0.25 dex metallicity bin in 
the interval -1.00 < [Fe/H] < +0.50, and a total of 20 dwarfs 
and subgiants of F5 to G8 spectral type (TableQ. 

2.2. Chiecks of sample contamination 

Our sample is restricted to objects closer than 40 pc. It is a very 
localised sample, c onsidering that the Galactic thin disk has a 
290 pc scale height dBuser et alJl998l) . This greatly reduces the 
probability of halo stars contaminating our sample, because the 
local density of halo stars is less than 0.05% (Buseret al. 1993). 
In order to reduce even further the probability of contamina- 
tion, we restricted our sample to relatively high metallicities 
([Fe/H] > -1.00). However, these criteria are not perfect, be- 
cause there exists an intersection between the metallicity distri- 
butions of disk and halo stars, i.e., there exists a non-negligible, 
albeit small, number of halo stars with [Fe/H] > - 1 .00. As an 
additional check of contamination, we performed a kinematic 
analysis of the objects. 

We developed a code to calculate the U, V, and W 
spatial velocity components . This code is based on the 
Johnson & Soderblom ('l987) formulation and makes use of 
the .Perrvman et al.. (.1998.) coordinate transformation matrix. 
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References: Coordinates: SIMBAD (FK5 system); parallaxes: Hipparcos catalogue i Perrvman & ESA 1997^; visual magnitudes: Bright Star 
Catalogue i Warren & Hoffleit 1987) for stars with an HR number and SIMBAD for those without; spectral types and lumin osity classes: 
Michigan Catalogue of HD Stars lliouk & Cowlevl[T97llHouklfT973.ll982tlHouk & Smith-Moordll988tlHouk & Swiftlll999l) for all stars, 
with the exception of HD 182 572 and HD 196 755 (Bright Star Catalogue), and HD 43 047 jFehrenbactll96ll) . 



It calculates the velocity components relative to the Sun, and 
then converts them to the local sta ndard of rest (LSR) using 
the velocity of the Sun taken from ' Pehnen & BinnevI (Il998l) . 
We adopted a right-handed coordinate system, in which U, V, 
and W are positive towards the Galactic center, the Galactic 
rotation and the Galactic north pole, respectively. The code in- 
put parameters were taken from the Hipparcos catalogue, with 
the exception of the radial velocities, which were determined 
by us when Doppler velocity corrections were applied to the 
FEROS spectra. Four of our objects (HD 22484, HD 22 879, 
HD 76 932, and HD 182 572) have ra dial velocities d etermined 
with the CORAVEL spectrometers dUdrv et al.lfl999i^ . which 
have high precision and accuracy, and there is excellent agree- 
ment between these and our results (linear correlation with dis- 
persion o- = 0.05 km s Mark that HD 182572, although 
initially selected, was not included in the final sample, be- 
cause we could not obtain Th spect ra for this object; nonethe- 
less, we used it for comparison with' Udrv et alJ ( 1 1999 *) because 
we had its FEROS spectrum and its derived radial velocity 
(-100.8 km s"'). Distances were calculated directly from the 
trigonometric parallaxes; extinction was not corrected for, as it 
can be considered negligible due to the proximity of our ob- 
jects (less than 40 pc). Our results are contained in Table |2 
along with the radial velocities used in the calculations. 

It is not possible to discriminate the disk from the halo 
population based on kinematic or metallicity criteria, indepen- 
dently. Nonetheless, this discri mination can be achi eved us- 
ing these criteria j'/wMZfaweoM.yZv. 'Schuster et'al] ( ll99 3^ demon- 
strated that halo stars whose metaUicities are higher than 



Table 2. Radial velocities (RV) and spatial velocity compo- 
nents (U, V, and W) of the sample stars, in a right-handed 
Galactic system and relative to the LSR. All values are in 

km s"^ 
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[Fe/H] = -1.00 have V velocity component lower than disk 
stars of the same metallicity, and that it is possible to sep- 
arate these populations with a line in the V vs. [Fe/H] di- 
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Fig. 1. V vs. [Fe/H] diagram for the sample stars. The diagonal 
line is a cut off between the halo and the disk populations (be- 
low and above, respectively). Velocity components V were cal- 
culated by us. The metallicities initially employed, represented 
as open squares, were taken from the literature (Table|3; small 
filled circles represent metallicities derived in this work. 



agram. Analyzing such diagram constructed for our sample 
stars (Fig. [Q, we can see that all of them are located far 
above the cut off line, indicating a highly remote probability 
of one of them belonging to the halo. Note that the metallic- 
ities used for the graph were taken from papers individually 
ch osen among the spectroscopic determinations that compose 
the lCavrel de Strobel et al ]|2001) catalogue, i.e., they were not 
determined by us. Even though they could have been wrong, 
these determinations would have to have been overestimated 
by at least 0.52 dex for a star to cross the line, which is very un- 
likely. After having determined our own metallicities (Table0, 
we confirmed that results from the literature are overestimated 
by 0.08 dex at most, when compared to our results. 

According to Buser et al] lll998h . the local density of thin 
disk, thick disk, and halo stars in the Galaxy conforms to the 
proportion 2000:108:1, respectively. If we take a completely 
random sample of 20 stars, there is a 34.9 % probability that it 
will be totally free of thick disk stars. The probability that it will 
contain 3 or more thick disk stars is low (8.0 %). Consequently, 
if there are thick disk stars in our sample, they can be regarded 
as a contamination without significant bearing on our results. 

2.3. ESO 1.52 m telescope spectra 

Observations with the ESO 1 .52 m telescope were performed 
by the authors in two runs, in March and August, 2001. All 
obtained FEROS spectra have high nominal resolving power 
(R = 48000), signal-to-noise ratio (S/N > 300 in the visible) 
and coverage (3500 A to 9200 A spread over 39 echelle or- 
ders). Spectra were acquired for all sample stars and for the 
Sun (Ganymede), and were used for atmospheric parameter and 
chemical abundance determination. Order #9 of the FEROS 
spectra, which extends from 4039 A to 4192 A, was used for 
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Fig. 2. Example of a solar (Ganymede) normalised FEROS 
spectrum (order #20), obtained with the ESO 1 .52 m telescope. 



the determination of Eu abundances by synthesis of the Eu n 
Une at 4129.72 A. 

FEROS spectra are reduced automatically by a MIDAS 
script during observation, immediately after the CCD read out. 
Reduction is carried out in the conventional way, applying 
the following steps: bias, scattered light, and flat field correc- 
tions, extraction, wavelength calibration (using ThAr calibra- 
tion spectra), and baricentric radial velocity correction. As a 
last step, the reduction script merges all echelle orders to form 
one continuous spectrum with a 5200 A coverage. However, the 
merge process is sometimes faulty, and strong curvatures and 
discontinuities are introduced in the final spectrum as result. 
It becomes difficult to distinguish between the "natural" cur- 
vature of the spectrum and that introduced by the merge, and 
the continuum normaUzation procedure is rendered very inse- 
cure. We chose to re-reduce all spectra, employing the same 
script used for the automatic reduction, but eliminating the 
merge, so that we could work with each order independently. 
Unfortunately, after the re-reduction we found out that each in- 
dividual order is plagued by a strong discontinuity itself, which 
demonstrates that some reduction steps prior to the merge are 
also faulty. We corrected this by simply trimming each order 
and retaining only the section beyond the discontinuity. As 
there are intersections in wavelength coverage between adja- 
cent orders, total loss of coverage caused by the trimmings was 
only ~ 500 A. After trimming the orders, we corrected them 
for the stellar radial velocity. Stars that were observed more 
than once had their spectra averaged, using their S/N ratio as 
weight. Finally, the spectra were normalised by fitting selected 
continuum windows with order 2 to 5 Legendre polynomials. 
Figure|2is an example of a solar (Ganymede) normalised ESO 
spectrum. 

2.4. OPD spectra 

Observations at OPD were performed by the authors in five 
runs, in May and October, 2000, and in May, August and 
October, 2002. The instrumental set-up was the following: 
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HD 



FEROS OPD CES + 3.60 m 



CES + CAT 



6500 6520 6540 6560 6580 6600 6620 

Fig. 3. Example of a solar (Ganymede) normalised OPD spec- 
trum. 

1800 lines mm"' grating, 250 )im slit, 1024 24)J.m-pixel back- 
illuminated CCD (#101 or #106). All obtained spectra have 
moderate nominal resolving power (R ~ 22000) and high 
signal-to-noise ratio (S/N > 200). The spectra are centered at 
the Ha line (6562.8 A) and cover ~150 A. Spectra were ob- 
tained for all sample objects, with the exception of HD 9562, 
and for the Sun (Ganymede), and were used for T^s determina- 
tion through profile fitting (Sect. ITT^ . FEROS spectra cannot 
be used for this purpose because the echelle order #29, which 
contains Ha, has a useful range of 6430 A to 6579 A, leaving 
only 16 A at the red wing of the line. Since Ha is very broad, 
it does not reach continuum levels this close to its center, pre- 
cluding normalisation. On the hand, OPD spectra reach 75 A at 
each side of Ha, allowing the selection of reliable continuum 
windows. 

Reduction was performed using the Image Reduction 
and Analysis Facility (IRAF^*), following the usual steps of 
bias, scattered light and flat field corrections, and extraction. 
Wavelength calibration was performed individually for each 
star by comparing the central wavelength of selected absorp- 
tion lines with their laboratory values. This way, Doppler cor- 
rections were not necessary, since the wavelength calibration 
corrects the spectra to the rest frame. Finally, spectra were av- 
eraged and normalised in the same way the ESO spectra were. 
Figure|3is an example of a solar (Ganymede) normalised OPD 
spectrum. 

2.5. ESO 3.60 m telescope and CAT spectra 

Observational spectra for the Th abundance determinations 
were obtained with the CES fed by the ESO 3.60 m telescope, 
for 16 of the 20 sample stars and for the Sun. Observations 
were carried out by the authors in January, 2002. Spectra were 
obtained with the blue optical path and the high resolution im- 



IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Europium Thorium 



2151 
9562 
16417 
20766 
20 807 
22484 
22 879 
30 562 
43 947 
52298 
59984 
63 077 
76932 
102 365 V 
128 620 V 
131117 V 
160 691 V 
196 378 V 
199 288 V 
203 608 V 



V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 



V 

V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 



V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 
V 



V 



V 
V 



V 



V 
V 
V 



V 



V 

V 
V 

V 



V 

V 
V 
V 
V 



age slicer (nominal resolving power R ~ 235 000). However, 
unknown problems during the observations degraded the effec- 
tive resolving power down to R ~ 130000. The spectra were 
centered at the Thii line (4019.13 A) and have a 27 A cover- 
age. Due to the low efficiency of the spectrograph at 4000 A, 
we limited the signal-to-noise ratio to S/N ~ 200 for stars with 
visual magnitude V > 6 and S/N ~ 300 for brighter objects. 

CES spectra were also obtained with the ESO's CAT, with 
the aim of determining the abundances of both Eu and Th. 
Observations were carried out by the authors in August, 1998. 
Not all sample objects have been observed: we obtained Eu 
spectra for 8 stars and the Sun, and Th spectra for 9 stars and 
the Sun. All spectra have high resolving power (R = 100000 
for Th, and R = 50 000 for Eu), and high signal-to-noise ratio 
(S/N ~ 300); coverage is 18 A. Table |3 details which objects 
were observed with which instruments and telescopes. 

Reduction followed the usual steps of bias and flat field 
corrections, extraction and wavelength calibration (using ThAr 
calibration spectra). Correction for the stellar radial velocities 
was implemented by comparing the wavelengths of a number 
of absorption lines with their rest frame wavelengths. Finally, 
the spectra were normalised by fitting selected continuum win- 
dows with order 2 to 5 Legendre polynomials. Stars that were 
observed in more than one night did not have spectra from dif- 
ferent nights averaged. Figure |4] shows as an example a solar 
(sky) normalised CES spectrum, observed with the 3.60 m. 
Examples of solar (sky) normalised CES spectra, observed with 
the CAT, are presented in Fig. [S] (for the Eu region), and in 
Fig m (for the Th region). 



3. Atmospheric parameters 

Atmospheric parameters (effective temperatures, surface grav- 
ities, microturbulence velocities (^), and metallicities) have 
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Fig. 4. Example of a solar (sky) normalised CES spectrum, ob- 
served with the ESO 3.60 m telescope, in the Thii line region 
at 4019. 13 A. The inset shows the Th ii line in greater detail. 



T3 




4122 4124 4126 4128 4130 4132 4134 4136 4138 
MA) 

Fig. 5. Example of a solar (sky) normalised CES spectrum, ob- 
served with the ESO CAT, in the Eun line region at 4129.72 A. 
The inset shows the Eu ii line in greater detail. 

been obtained through an iterative and totally self-consistent 
procedure, composed of the following steps: 

1. Determine photometric T^s ([Fe/H] from literature). 

2. Determine log ^. 

3. Determine Teff from Ha fitting. 

4. Set Teif = (photometric + Hq')/2. 

5. Determine [Fe/H] and ^. 

6. Determine photometric Teff ([Fe/H] from step 5). 

7. Set Teff = (photometric + Hq')/2. 

8. Have the parameters changed since step 8 of last iteration? 
If yes, go back to step 2. 

9. END. 

Notice that a minimum of 2 iterations is required, since in 
step 8 we must compare the results obtained with those from 
the preceding iteration. Input parameters for each step are sim- 
ply the output parameters of the preceding step. At all times. 



4010 4012 4014 4016 4018 4020 4022 4024 4026 
X(A) 

Fig. 6. Example of a solar (sky) normalised CES spectrum, ob- 
served with the ESO CAT, in the Th n line region at 4019. 13 A. 
The inset shows the Th ii line in greater detail. 



the adopted effective temperature is the average of the two de- 
terminations (photometric and from Ho- fitting). Owing to this, 
a new average is performed each time a new determination is 
made (steps 4 and 7, following Tetf determinations by Wa fit- 
ting and photometry, respectively). The techniques used at each 
step are described in the next sections, and the final results are 
presented in Table 



3A. Model atmospheres and partition functions 

The determination of effective temperatures through Ha profile 
fitting, and of microturbulence velocities and chemical abun- 
dances through detailed spectroscopic analysis, demands the 
use of model atmospheres. Those employed by us were the 
plane parallel, flux constan t, LTE model atmosphere s described 
and discussed in detail bv'E dvardsson et alJ ( Il993l) . The T(t) 
laws were interpolated using a code gently made available 
by Monique Spite (Observatoire de Paris-Meudon, France). 
Throughout this work, we have used [Fe/H] as the metallicity 
of the model atmosphere. Partition functions were calculated 
by us using the polynomial approximations of llrwinl ( ll98ll) . 



3.2. Effective temperature 
3.2.1 . Photometric calibrations 



We determined effective temperatures using the photomet- 
ric calibrations from P orto de Mello et alJ (|2005). These 
metallicity-dependent calibrations were constructed using 
highly accurate effective temperatures obtained with the in- 
frared flux method (iBlackwell et al..l98Q..1991.) . and make use 
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of the (B - V) and (V — K) Johnson, ib - y) and p Stromgren, 
and (Bt-Vj) Hipparcos colour indices: 

reff(K) = 7747 - 3016 (B - V) {1 - 0.15 [Fe/H]), cr = 65 K 

Teff (K) = 8974 - 2880 (V - K) + 440 (V - K)2, cr = 50 K 

reff(K) = 8481 -6516 (b-y){l -0.09 [Fe/H]), cr = 55K 

reff(K) = 1 1 654 Vj6 - 2.349, cr = 70 K 

reff(K) = 7551 - 2406 (Bt - Vt) {1 - 0.2 [Fe/H]), cr = 64 K 

A comple te description of the derivat ion of these calibrations 
IS given in 'P orto de Mello et alJ ( Eool . 

Table 0] contains photometric data from the literature. 
The only sample star absent on the table is HD 128 620 
(a Cen A), whose photometric data is unreliable due to its 
high apparent magnitude. The metallicities presented in the 
table were taken frorn paper s individually chosen from the 
ICavrel de Stro bel et al.' ('200 1') catalogue, with the exception of 
HD 160691 (Porto de Mello et al. 2005), and were only used 
on step 1 of the atmospheric parameter det ermination for each 
star. The (b - y) colour indices taken from lGr0nbech & OlsenI 
( 1976, 1977) and Crawford et al. ( 197o|) were originally in the 
loiseni (119831) photomet ric sys t em. S ince the calibration used 
by us was based on the lOlsenI ( Il993l) photometric system, we 
converted the afore-mentioned data to this system through the 
relation (fe-y)oisen i92i^0.8858(^-y)oisen 1983+0.0532, which 
is Equation (1) from^O lsenl ( Il993l) . 

Photometric temperatures have been averaged using the 
maximum likelihood method, assuming that each value follows 
a Gaussian probability distribution, which results in 




where is the number of calibrations effectively used. We took 
the standard deviation of the averages as their uncertainties. 
These vary from 27 K to 37 K, and such a short range prompted 
us to adopt a single value cr = 32 K for the average photometric 
Teff of each sample star Take heed that th ese are internal uncer- 
taintie s only; for a detailed discussion, see lPorto de Mello et alJ 
(12005'). Final photometric effective temperatures are presented 
in TableEl 



Table 5. Uncertainty of the effective temperatures determined 
by He profile fitting, using the Sun as a standard. For detailed 
explanation, see text. 



Parameter 


Parameter 


Tetf 




uncertainty 


uncertainty 


Continuum 


0.20 % 


25 K 




0.20 dex 


20 K 


[Fe/H] 


0.08 dex 


5K 




0.10 km S-' 


5K 


Personal judgement 




20 K 


Blue/red wings 




20 K 


TOTAL 




43 K 



In order to check the accuracy of the determinations, we 
analysed the solar spectrum (Fig. |7J. The adopted solar pa- 
rameters were T^g = 5777 K (Neckel 1986), \ogg = 4.44, 
[Fe/H] - 0.00, rnicrot urbulence velocity ^ - 1.15 km s ' 
(lEdvardsson eTaP Il993h . and n(He)/n(H) = 0.10. The ob- 
tained temperature was 5767 K, 10 K lower than the solar value 
adopted by us. This small discrepancy was corrected by adding 
-H 10 K to all stellar effective temperatures determined by Ha 
profile fitting. 

Uncertainties in the determinations were estimated using 
the Sun as a standard (Table |5}. The influence of atmospheric 
parameter uncertainties were evaluated changing these param- 
eters (one at a time) by an amount deemed representative of 
their uncertainties, and then redetermining the effective tem- 
perature. We tested the influence of the continuum placement 
uncertainty, which is ~0.2%, by multiplying the flux by 1.002, 
and redetermining the temperature. Since the fitting of theoret- 
ical profile to the observation is done via a. fit-by-eye technique, 
the personal judgement of who is making the analysis intro- 
duces uncertainty to the obtained temperature; we estimated 
this uncertainty to be ~20 K. Finally, some uncertainty comes 
from the difference between blue and red wings, which are best 
fit by profiles of effective temperatures differing by ~20 K. 
The total uncertainty was estimated at 43 K by the root-mean- 
square (RMS) of the individual sources cited above. Final Ha 
effective temperatures are presented in Table 



3.2.2. Ha profile fitting 

The Hg profiles have been studied in deta il (bvlda Silva|l975| 
lGehr en'l98l'. 'Fuhrmann et al."l993', 'l994', iGratton e t al.'l996'. 
and ICowlev & Castelli 2002i among others), and have been 
shown, in the case of cool stars, to be rather insensitive to sur- 
face gravity, microturbulence velocity, metallicity, and mixing 
length parameter variations. They are, notwithstanding, very 
sensitive to the effective temperature of the atmosphere. By 
comparing the observations with theoretical profiles calculated 
for various effective temperatures, we can estimate this atmo- 
spheric parameter. 

Theoretical profiles were calculated with a code devel- 
oped by us from the original routines of Pr aderie (1_967)- 
This code takes into account the convolution of radia- 
tive. Sta rk (Vidal et al." '1971*), Doppler and self-resonance 
JCavrel & Traving.l96Q,) broadenings. 



3.2.3. A(dopted effective temperatures 

The adopted affective temperature for each star was the arith- 
metic mean of the photometric and Ha determinations, after 
convergence was achieved in the iterative procedure. The un- 
certainty of the adopted Teff was taken as the standard deviation 
of the mean, resulting in 27 K. 

3.3. Surface gravity 

Suiface gravities were estimated from effective temperatures, 
stellar masses and luminosities, employing the known equation 

log(^) = log(— ) +41og(|^) + 0.4(Mboi-Mboio), (1) 

in which we adopted Mboio - 4.75 (lNeckeJl986l) . 
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Table 4. Photometric indices and metallicities for all sample stars. Metallicities were taken from the literature for use in the 
kinematic characterization of the sample (Sect. 12. 2> and as input for the first step in the iterative determination of atmospheric 
parameters (Sect. 13. 2t . 



HD 


(B-V) 


Ref. 


(b-y) 


Ref. 


(V-K) 


Ref. j8 


Ref. 


(Bt-Vt) 


Ref. 


[Fe/H] 


Ref. 


2151 


0.62 


1 


0.379 


3 


— 


2.597 


11 


— 


— 


-0.11 


15 


9562 


0.64 


1 


0.408 


3 


1.422 


13 2.585 


10 


0.709 


14 


+0.13 


16 


16417 


0.66 


1 


0.412 


3 


— 


— — 


— 


0.730 


14 


+0.00 


17 


20766 


0.64 


1 


0.404 


4 


1.537 


13 2.586 


5 


— 


— 


-0.22 


18 


20 807 


0.60 


1 


0.383 


4 


— 


2.592 


5 


— 


— 


-0.22 


18 


22484 


0.58 


1 


0.363 


3 


1.363 


13 2.608 


5 


0.626 


14 


-0.13 


19 


22 879 


0.54 


2 


0.369 


3 


— 


— — 


— 


0.581 


14 


-0.84 


20 


30562 


0.62 


1 


0.403 


3 


1.410 


13 2.610 


5 


0.709 


14 


+0.27 


21 


43 947 


— 


— 


0.377 


6 


— 


2.598 


6 


0.604 


14 


-0.33 


19 


52298 


0.46 


2 


0.320 


3 


— 


— — 


— 


0.500 


14 


-0.84 


22 


59 984 


0.54 


1 


0.354 


7 


— 


2.599 


5 


0.566 


14 


-0.75 


23 


63 077 


0.60 


1 


0.387 


5 




2.590 


5 






-0.78 


23 


76932 


0.53 


1 


0.368 


5 


1.410 


13 2.595 


5 


0.556 


14 


-0.76 


16 


102 365 


0.66 


1 


0.411 


3 




2.588 


6 






-0.28 


24 


131 117 


0.60 


1 


0.389 


4 




2.621 


9 


0.662 


14 


+0.09 


16 


160 691 


0.70 


1 


0.433 


3 








0.786 


14 


+0.27 


24 


196 378 


0.53 


1 


0.369 


5 




2.609 


5 


0.579 


14 


-0.44 


23 


199 288 


0.59 


2 


0.393 


6 




2.588 


6 


0.638 


14 


-0.66 


25 


203 608 


0.49 


1 


0.338 


8 


1.310 


12 2.611 


8 


0.522 


14 


-0.67 


23 


References: 
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h 2 


- j^ermilliod 


1 1987); 3 - Olsen ( 1994); 4 - Olsen 


(1991; 


5 - lOrOnbech & 



6 - Olsen 1 1983); 7 - Manfroid & Sterken ( 1987); 8 



11 - i Hauck & Mermilliod 1 1998); 12 - Koornneei 41983 
Jl999l): 16 -iGr attonetal. 1 1996); 17 - Gehre^ 



13 - Idi Benedettd jl998h : 1 4 - .Perrvman & ESAi |1997^: 15 - ICastro et all 



da Silva & Por to de Melldl 1200fl); 22 - IHartmann & GehrerJ 
Axeretalj ^99^ 



23 - lEdvardsson et alj 1199J; 24 - iPorto de Mello et al 



Idel Peloso et alj llOOd): 19 - IChen et alJ J200(j): 20 - iFuhrmanrJ jl99^: 21 




MA) 

Fig. 7. Ha profile of the Sun. Dashed line is the best fit 
(5767 K). Dotted lines show profiles computed when the ef- 
fective temperature was changed by +50 K and +100 K. 



Luminosities were calculated by us from the visual 
magnitudes of T able \l\ and bolometric corrections from 
iHabets & Heintzeitl981 ). Stellar masses were estimated using 
HR diagrams with the Gen92/96 sets of evolutionary tracks. 
Each diagram was calculated by Gen92/96 for 1 of 5 different 
metallicities (Z = 0.0010, 0.0040, 0.0080, 0.0200, and 0.0400, 
with Zq - 0.0188), with tracks for 5 diff'erent masses (0.80, 
0.90, 1.00, 1.25, and 1.50 m©). Masses were estimated for each 
star, interpolating among the tracks calculated for the metallici- 




Fig. 8. Example of a Gen92/96 set of evolutionary tracks, for 
Z = 0.0200. Filled squares are the sample stars, after conver- 
gence of the atmospheric parameter determination procedure. 
The Sun is represented by the symbol O. The effective temper- 
atures and metallicities are the final, adopted values (Table0. 



ties above. Then, we interpolated among metallicities to obtain 
the correct mass for the star Figure |8l presents one such dia- 
gram, calculated for Z - 0.0200, in which we used the final, 
adopted effective temperatures and metallicities (TableC). 



E.F. del Peloso et al.: The age of the Galactic thin disk from Th/Eu nucleocosmochronology. I 



9 



In order to check the accuracy of the determinations, we 
estimated the solar mass itself in the Gen92/96 diagrams. The 
obtained value was 0.97 ihq. This discrepancy was corrected by 
adding +0.03 wjq to the masses obtained for the sample stars. 

The uncertainty of the surface gravities was obtained by 
error propagation from Equation Q, resulting in the expression 



O-log 



ig(m/m0) 



16 (T? ^ ; 

log r,.(r' 



analysis showed that the influence of the bolometric magnitude 
uncertainty is negligible, and as such it was not taken into con- 
sideration. 

Stellar mass uncertainty is a function of the position of the 
star in the HR diagram. Stars located in regions where the evo- 
lutionary tracks are closer together will have larger errors than 
those located between tracks further apart. Nevertheless, we 
have chosen to adopt a single average uncertainty of 0.03 m©, 
so as to simplify the calculations. Taking stellar masses and 
effective temperatures typical of our sample (1.00 and 
6000 K, respectively), we have ixiog 7-^^ = 1.9 lO^^* and 
o"iog(m/mQ) - 1.3 10"^, and consequently cr\og g - 0.02 dex. 
This value does not take into consideration the uncertainties in- 
trinsic to the calculation of evolutionary tracks, which are very 
diflicult to assess. 

3.4. Microturbulence velocity and metallicity 

Microturbulence velocities and metallicities have been deter- 
mined through detailed, diff'erential spectroscopic analysis, rel- 
ative to the Sun, using EWs of Fe i and Fe n lines. In such differ- 
ential analysis, the Sun is treated like any other sample star: its 
spectra are obtained under the same conditions, with the same 
resolving power, signal-to-noise ratio, and equipment, and are 
reduced by the same procedures; the same model atmospheres 
and analysis codes are utilised. This way, systematic errors and 
eventual non-LTE effects are reduced, because they are par- 
tially canceUed out. 

3.4.1. Equivalent width measurements 

We were not interested solely in Fe lines, but also in the lines 
of other elements that contaminate the spectral regions of the 
Thii line in 4019.13 A (V, Cr, Mn, Co, Ni, Ce, and Nd), and 
of the Euii line in 4129.72 A (Ti, V, Cr, Co, Ce, Nd, and Sm), 
whose EWs were then used for the abundance determinations 
of Sect, a 

An initial list of l i nes was constructed from iMoore et alJ 
("1966'), 'Steffen' (1985'), Cavrel de Strobel & Bentolila ( 1989), 
Brown & Wallerstein ( .1992) , Furenlid & Mevlan (1990), and 
Mevlan et alJ (Il993ir Then, a series of selection criteria were 
applied: lines with wavelengths lower than 4000 A were re- 
moved, since the line density in this region is very high, low- 
ering the continuum and rendering its normalization insecure, 
besides making it too difficult to find acceptably isolated lines; 
lines with wavelengths greater than 7000 A were discarded, 
due to excessive telluric line contamination in this region; lines 
with EW > 1 10 mA in the Sun, which are more sensitive to un- 
certainties in microturbulence velocity due to saturation, were 
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Fig. 9. Panel (a): Compari son between our solar EWs and those 
from 'Mevlan e talJ(ll993h . Solid line corresponds to a y - x 
relation, and the dotted one is a linear regression. Panel (b): 
Percent differences between the same EWs. Solid lines corre- 
spond to mean differences, and the dotted ones represent the 
mean differences modified by their standard deviations (±1 cr). 



eliminated; only isolated lines were kept, eliminating those 
with known contaminations; lines whose EWs could not have 
been, for any reason, accurately measured in the solar spectra, 
were rejected. 

EWs have been measured by Gaussian profile fitting. In or- 
der to check the measure ment accuracy, we c ompared our so- 
lar EWs with fliose from lMevlan et alJ ( Il993l) . The EWs from 
'Mevlan et al." were measured by Voigt profile fitting, using the 



JCur ucz et al. ( 1984) solar atlas, which has very high resolving 
power and signal-to-noise ratio (R = 522 000 and S /N - 3000, 
between 4500 A and 6450 A). A linear fit between the two data 
sets. 



EW 



Meylan et al, 



1993 - (1.03463 + 0.00649) EWthis work? 



shown in panel (a) of Fig. |5| presents low dispersion (cr = 
2.9 mA) and high correlation coefficient {R - 0.994). This re- 
lation is used to correct all our EWs. The mean percent differ- 
ence between the measurements is -2.5%, presenting no de- 
pendence on the EW (Fig. |9] panel (b)). The larger percent 
differences exhibited by the weakest lines are expected, since 
these are more prone to uncertainties. More than 86% of the 
lines agree at a 10% level. 
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Table 6. A sample of the EW data. The complete content of 
this table is only available in electronic form at the CDS. For a 
description of the columns, see text (Sect. imi . 



A (A) 


Element 


;r(eV) 




Sun ■■ 


• HD 203 608 


5657.436 


Vi 


1.06 


-0.883 


9.5 ■■ 


0.0 


5668.362 


Vi 


1.08 


-0.920 


8.5 •■ 


0.0 


5670.851 


Vi 


1.08 


-0.452 


21.6 •■ 


0.0 


5427.826 


Fell 


6.72 


-1.371 


64 •■ 


0.0 


6149.249 


Fell 


3.89 


-2.711 


40.9 • ■ 


24.7 



A last selection criterium was applied to the line list: given 
the 2.9 mA dispersion of the l inear relation between our EWs 
and those from M evlan et all we rejected all Fe lines with 
EW < 6 mA (i.e., ~ 2cr). This criterium was not implemented 
for the lines of the other elements because they have fewer in- 
teresting lines available, and a cut off at 6 mA would be too re- 
strictive. The final list of lines is composed of 65 Fe i, 10 Fe ii, 
31 Ti, 6 V, 21 Cr, 6 Mn, 8 Co, 13 Ni, 5 Ce, 2 Nd, and 1 Sm Hnes, 
in the Sun. Not all these lines have been used for all sample 
stars, since random events, like excessive profile deformation 
by noise, prevented sometimes an accurate measurement. 

Table|6lpresents a sample of the EW data. Its complete con- 
tent, composed of the EWs of all measured lines, for the Sun 
and all sample stars, is only available in electronic form at the 
CDS."* Column 1 lists the central wavelength (in angstroms). 
Column 2 gives the element symbol and degree of ionization. 
Column 3 gives the excitation potential of the lower level of 
the electronic transition (in eV), Column 4 presents the solar 
log gf derived by us, and the subsequent Columns present the 
EWs, in mA, for the Sun and the other stars, from HD 2151 to 
HD 203 608 (in order of increasing HD number). 

3.4.2. Spectroscopic analysis 

The first step of the spectroscopic analysis consisted of deter- 
mining solar log gf values for all measured lines. For such 
determination, we employed a code written by us, based on 
routines originally developed by Monique Spite. This code 
calculates the \og gf that an absorption line must have for 
its abundance to be equa l to the standard solar value from 
iGrevesse & Sauvail ( Il998h . given the solar EW and lower ex- 
citation potential of the line. For all stars, chemical abundances 
were determined using these solar log gf values. 

Microturbulence velocity and metallicity have been deter- 
mined through a procedure which is iterative, and which can be 
summarised by the following sequence of steps; 

1 . Determine [Fe/H] for all lines, using 

[Fe/H] model atmosphere ~ [Fe/H] literature- 

2. Set [Fe/H]modei atmosphere = average [Fe/H] from lines. 

3. Determine [Fe/H] for all lines. 

4. Is [Fe/H]modei atmosphere = average [Fe/H] from lines? 
If no, go back to step 2. 

* Via anonymous ftp to cdsarc.u-strasbg.fr or via |http://| 
cdsweb .u-strasbg . fr/cgi-bin/qcat? 



5. Draw a [Fe/H] vs. EW graph for all lines, and fit a line. 

6. Is the angular coefficient of the line negligible? 
If no, modify ^ and go back to step 3. 

7. END. 

In steps 1 and 3, iron abundances are calculated line-by- 
line using a code very similar to the one used for the determi- 
nation of the solar log gf values. The only difference is that, 
in this case, the log gf are fixed at their solar values, and an 
abundance is calculated for each line. In step 6, the angular co- 
efficient is considered negligible when it is less than 25% of 
its own uncertainty (obtained by the fitting calculation). It is 
important to note that effective temperature and surface gravity 
remain constant at all times. In practice, the procedure is exe - 
cuted automatically by a code, developed by 'da Silva' 1200(f). 
This code iterates the parameters, calculating abundances and 
linear fits, and making the appropriate comparisons, without 
human interference. 

Iron abundances obtained using Fei and Fen lines agree 
very well, their differences being always lower than dispersion 
among the lines (as can be seen in Table |8}. The differences 
show no dependence on any atmospheric parameter Ionization 
equilibrium was not enforced by the analysis procedure, but it 
was, nonetheless, achieved. 

Microturbulence velocity uncertainties have been estimated 
by varying this parameter until the uncertainty of the angular 
coefficient of the linear fit to the [Fe/H] vs. EW graph became 
equal to the coefficient itself. We found out that the uncertain- 
ties are metallicity-dependent, and can be divided into three 
ranges: 

("0.23 km s"', if [Fe/H] < -0.7 

0-^ = 1 0.13kms"', if -0.7 < [Fe/H] < -0.3 
[ 0.05 km s ', if -0.3 < [Fe/H] 

This behaviour is expected, as the dispersion between Fe lines 
increases for metal-poor stars. Metallicity uncertainties were 
estimated in the same way that uncertainties for the other ele- 
ments were (Sect. 14.21 . and are presented in Table0 along with 
the final, adopted values of all atmospheric parameters. 

4. Abundances of contaminating elements 

4.1. Abundance determination 

The abundance of the elements other than iron have been deter- 
mined using the same code as in Sect. 13.4.21 with the adopted 
atmospheric parameters and the solar log gf values. 

Adopted abundances have been obtained taking the arith- 
metic mean of the abundances from all lines of each element, 
irrespective of ionization state. When dispersion between lines 
was less than 0.10 dex, we eliminated those that deviated by 
more than 2cr from the mean, and when it was more than 
0.10 dex, we eliminated those that deviated by more than 1 cr. 
The number of lines kept is presented in Tables |8l and |9] along 
with adopted abundances, dispersions between the lines, and 
uncertainties, calculated according to the procedure described 
in Sect. lO 
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Table 7. Adopted atmospheric parameters, including the pho- 
tometric and Ha effective temperatures used to obtain the 
adopted mean values, and the stellar masses used to obtain the 
surface gravities. Uncertainties of photometric. Ha, and mean 
effective temperatures, and of stellar masses and surface gravi- 
ties, are the same for all stars: 32 K, 43 K, 27 K, 0.03 m©, and 
0.02 dex, respectively. 



HD 




Ta (K) 


mo 




[Fe/H] 


f (km s"') 




Phot. 


Ha MEAN 














2151 


5909 


5799 


5854 


1.19 


3.98 


-0.03 


±0.09 


1.32 


±0.05 


9562 


5794 


— 


5794 


1.23 


3.95 


+0.16 


±0.09 


1.45 


±0.05 


16417 


5821 


5817 


5819 


1.15 


4.07 


+0.13 


±0.09 


1.38 


±0.05 


20766 


5696 


5715 


5706 


0.95 


4.50 


-0.21 


±0.09 


1.01 


±0.05 


20 807 


5866 


5863 


5865 


0.99 


4.48 


-0.23 


±0.10 


1.18 


±0.05 


22484 


5983 


6063 


6023 


1.16 


4.11 


-0.03 


±0.09 


1.44 ±0.05 


22 879 


5928 


5770 


5849 


0.73 


4.34 


-0.76 


±0.12 


0.69 


±0.23 


30562 


5883 


5855 


5869 


1.18 


4.09 


+0.19 


±0.09 


1.51 


±0.05 


43 947 


5940 


5937 


5889 


0.95 


4.32 


-0.27 


±0.10 


1.02 


±0.05 


52298 


6305 


6200 


6253 


1.09 


4.41 


-0.31 


±0.10 


1.44 


±0.13 


59 984 


5968 


5848 


5908 


0.97 


3.96 


-0.67 


±0.11 


1.07 


±0.13 


63 077 


5752 


5713 


5733 


0.77 


4.15 


-0.76 


±0.11 


0.78 


±0.23 


76932 


5874 


5825 


5850 


0.82 


4.14 


-0.84 


±0.11 


0.84 


±0.23 


102 365 


5705 


5620 


5663 


0.88 


4.43 


-0.29 


±0.09 


1.05 


±0.05 


128 620 




5813 


5813 


1.09 


4.30 


+0.26 


±0.09 


1.23 


±0.05 


131 117 


5994 


5813 


5904 


1.25 


3.96 


+0.10 


±0.09 


1.49 


±0.05 


160 691 


5736 


5676 


5706 


1.06 


4.19 


+0.28 


±0.09 


1.27 


±0.05 


196 378 


6014 


6044 


6029 


1.14 


3.97 


-0.37 


±0.10 


1.64 


±0.13 


199 288 


5785 


5734 


5760 


0.82 


4.35 


-0.59 


±0.10 


0.84 


±0.13 


203 608 


6057 


5987 


6022 


0.86 


4.31 


-0.67 


±0.11 


1.18 


±0.23 


The 


lines 


of V, 


Mn, 


and Co exhibit hyperfine 


struc- 



tures (HFSs) strong enough that they must be taken into ac- 
count, or else severe inaccuracies would be introduced. At 
the resolving power achieved by FEROS, line profiles re- 
main Gaussian even for elements with important HFSs, so 
that we can still use the EWs we measured by Gaussian pro- 
file fitting. For V and Mn lines, we took the wavelengths o f 
the HFS components from the line Usts of iKuruc^ ( l2003l) . 
whereas for t he Co lines, we c alculated the wavelengths our- 
selves ( del Peloso et al.lEo05t ). These calcu lations were real- 



ized using Cas imir's equation i CasirniJ 19631) and HFS labora- 
tory data from'Guthohrlein & Keller' (" 1990*), 'Pickering ('l996'), 
and Picker ing & Thome (1996J. The relative intensities of the 
components were ob tained from the calculatiori s of White & 
Eliason, tabulated in ICondo n & Sho rtlevI ill96 7^. A modified 
version of the code, which takes HFS into consideration, was 
used for the determination of the log gf values and abundances. 



4.2. Uncertainty assessment 

Chemical abundances determined using EWs are subject to 
three sources of uncertainty: atmospheric parameters, EWs, 
and model atmospheres. If these sources were independent 
from one another, the total uncertainty could be estimated by 
their RMS. However, there is no such independency, as can be 
clearly perceived when the iterative character of the analysis is 
called to mind. In this circumstance, RMS yields an estimate 
of the maximum uncertainty. Given the extreme complexity of 



6300 



6100 - 



6000 



5900 



5800 - 



HD 22484 



HD 63077 



HD 160691- 



-0.8 -0.6 -0.4 -0.2 0.0 0.2 

[Fe/H] 

Fig. 10. Stars used as uncertainty standards are labelled. 
HD 203 608 was not chosen when assessing the uncertainties of 
contaminating elements, because it lacks V and Nd abundance 
determinations; however, it was chosen as the hot, metal-poor 
standard when assessing Eu and Th uncertainties. The effec- 
tive temperatures and metallicities are the adopted ones. Dotted 
lines are an example of the weights used for uncertainty deter- 
mination in stars not chosen as standards - see text for details. 



the interdependency between the sources of uncertainty, we did 
not attempt to disentangle them, and used the RMS instead. 

We chose four stars as standards for the assessment, with 
extreme effective temperatures and metallicities: HD 160691 
(cool and metal-rich), HD 22484 (hot and metal-rich), 
HD 63 077 (cool and metal-poor), and HD 59984 (hot and 
metal-poor). HD 203 608 is as metal-poor as, but hotter than 
HD 59 984, and would apparently be better suited as a standard 
star; yet we did not use it, because it lacks measured V and Nd 
lines, and thus could not be used to estimate the uncertainties 
for these elements (Fig.llO>. 

Initially, we determined the iron abundance uncertainty. 
For this purpose, we evaluated the influence of all uncertainty 
sources, independently. The influence of effective temperature, 
surface gravity and microturbulence velocity was obtained by 
recalculating the Fe abundances using the adopted atmospheric 
parameters modified, one at a time, by their own uncertainties: 
+27 K, +0.02 dex, and +0.05 km s"' (metal-rich stars) and 
+0.23 km s"' (metal-poor stars), respectively. The difference 
between the adopted Fe abundance and the one obtained with 
the modified parameter is our estimate of the influence of the 
parameter on the total uncertainty. 

The influence of EWs was determined by recalculating the 
abundances using EWs modified by percentages obtained when 
comparing our measurements with those from 'Mevlan e t alJ 
( 1993) (Fig.|9}. This comparison allows us to estimate the influ- 
ence of the EW measurement per se, as well as the influence of 
continuum placement. The percent difference between the two 
sets of measurements presents a 15% dispersion for lines with 
EW < 20 mA, and 5% for those with EW > 20 mA. The Fe 
EWs were, then, modified by these amounts. However, this was 
only done for the two metal-rich standard stars. For the metal- 



Table 8. Fe, Ti, V, Cr, and Mn abundances, relative to H. N is the number of absorption lines effectively used for each abundance determination, and ct^x/h] is the 
deviation of the lines of element X. The uncertainties of abundances relative to H and to Fe, calculated as described in Sect. 14.21 are also presented (as uncert. 
uncert.[x/Fe], respectively). 



HD 


[Fei/H] 


N 


cr [Fel/Hl 


[Feii/H] 


N 


O" [Fcll/Hl 


[Fe/H] 


N 


C [Fe/H] 


uncert.iFc/H] 


[Ti/H] 


N 


C [Ti/Hl 


uncert.iTi/Hj 


uncert.[Ti/Fe] 


2151 


-0.03 


54 


0.08 


-0.03 


8 


0.05 


-0.03 


62 


0.08 


0.09 


+0.00 


29 


0.12 


0.09 


0.02 


9562 


+0.16 


55 


0.06 


+0.19 


9 


0.04 


+0.16 


64 


0.06 


0.09 


+0.15 


29 


0.09 


0.09 


0.02 


16417 


+0.13 


55 


0.04 


+0.12 


8 


0.03 


+0.13 


63 


0.04 


0.09 


+0.14 


28 


0.06 


0.09 


0.02 


20766 


-0.21 


51 


0.05 


-0.21 


8 


0.05 


-0.21 


59 


0.05 


0.09 


-0.16 


29 


0.06 


0.09 


0.02 


20 807 


-0.23 


48 


0.04 


-0.25 


8 


0.04 


-0.23 


56 


0.04 


0.10 


-0.18 


28 


0.07 


0.09 


0.02 


22484 


-0.03 


51 


0.07 


-0.03 


8 


0.05 


-0.03 


59 


0.07 


0.09 


+0.01 


28 


0.09 


0.10 


0.01 


22 879 


-0.75 


38 


0.12 


-0.80 


8 


0.08 


-0.76 


46 


0.12 


0.10 


-0.46 


27 


0.10 


0.09 


0.02 


30 562 


+0.19 


58 


0.05 


+0.19 


7 


0.04 


+0.19 


65 


0.05 


0.09 


+0.22 


28 


0.06 


0.09 


0.02 


43 947 


-0.28 


49 


0.06 


-0.22 


8 


0.06 


-0.27 


57 


0.07 


0.10 


-0.23 


27 


0.10 


0.09 


0.02 


52 298 


-0.31 


38 


0.06 


-0.28 


8 


0.06 


-0.31 


46 


0.06 


0.10 


-0.23 


25 


0.12 


0.09 


0.02 


59 984 


-0.67 


43 


0.10 


-0.68 


8 


0.11 


-0.67 


51 


0.01 


0.11 


-0.48 


28 


0.12 


0.08 


0.01 


63 077 


-0.76 


43 


0.12 


-0.78 


8 


0.10 


-0.76 


51 


0.11 


0.09 


-0.48 


25 


0.09 


0.09 


0.03 


76932 


-0.84 


36 


0.11 


-0.84 


8 


0.10 


-0.84 


44 


0.11 


0.10 


-0.57 


27 


0.08 


0.09 


0.02 


102 365 


-0.28 


52 


0.05 


-0.33 


8 


0.05 


-0.29 


60 


0.05 


0.09 


-0.17 


28 


0.05 


0.09 


0.02 


128 620 


+0.26 


56 


0.06 


+0.26 


9 


0.04 


+0.26 


65 


0.06 


0.09 


+0.29 


27 


0.08 


0.09 


0.02 


131 117 


+0.09 


53 


0.05 


+0.14 


9 


0.05 


+0.10 


62 


0.05 


0.09 


+0.09 


29 


0.06 


0.09 


0.02 


160691 


+0.27 


52 


0.05 


+0.33 


8 


0.03 


+0.28 


60 


0.05 


0.09 


+0.30 


29 


0.07 


0.09 


0.02 


196 378 


-0.37 


32 


0.08 


-0.33 


9 


0.06 


-0.37 


41 


0.07 


0.10 


-0.30 


20 


0.08 


0.09 


0.01 


199 288 


-0.59 


37 


0.07 


-0.59 


8 


0.06 


-0.59 


45 


0.06 


0.10 


-0.38 


22 


0.09 


0.09 


0.02 


203 608 


-0.68 


36 


0.10 


-0.60 


9 


0.11 


-0.67 


45 


0.11 


0.10 


-0.62 


20 


0.12 


0.09 


0.02 



HD 


[V/H] 


N 


cr [v/H] 


uncert. [V/H] 


uncert. [v/Fe] 


[Cr/H] 


N 


o- [aim 


uncert. [c/H] 


uncert. [Cr/Fe] 


[Mn/H] 


N 


C [Mn/HJ 


uncert. [Mn/H] 


uncert.[Mn/Fc] 


2151 


+0.02 


4 


0.01 


0.10 


0.04 


-0.04 


19 


0.06 


0.10 


0.02 


-0.09 


6 


0.04 


0.08 


0.03 


9562 


+0.16 


6 


0.07 


0.11 


0.04 


+0.18 


20 


0.07 


0.10 


0.02 


+0.16 


6 


0.05 


0.07 


0.02 


16417 


+0.18 


6 


0.07 


0.11 


0.04 


+0.14 


19 


0.04 


0.10 


0.02 


+0.14 


6 


0.03 


0.07 


0.03 


20 766 


-0.25 


5 


0.07 


0.11 


0.04 


-0.21 


19 


0.06 


0.10 


0.02 


-0.30 


5 


0.03 


0.08 


0.03 


20 807 


-0.23 


4 


0.03 


0.11 


0.04 


-0.23 


19 


0.08 


0.11 


0.02 


-0.29 


6 


0.03 


0.08 


0.03 


22484 


+0.06 


4 


0.08 


0.06 


0.03 


-0.03 


20 


0.07 


0.10 


0.02 


-0.07 


6 


0.03 


0.09 


0.04 


22 879 


-0.50 


4 


0.15 


0.13 


0.05 


-0.69 


13 


0.10 


0.11 


0.01 


-0.96 


3 


0.10 


0.08 


0.03 


30 562 


+0.23 


6 


0.10 


0.10 


0.04 


+0.20 


19 


0.08 


0.10 


0.02 


+0.18 


6 


0.02 


0.07 


0.03 


43 947 


-0.12 


4 


0.08 


0.11 


0.04 


-0.27 


20 


0.08 


0.11 


0.02 


-041 


6 


0.09 


0.08 


0.03 


52 298 


-0.26 


2 


0.09 


0.10 


0.04 


-0.30 


17 


0.13 


0.11 


0.02 


-0.39 


6 


0.09 


0.08 


0.03 


59 984 


-0.44 


4 


0.08 


0.13 


0.06 


-0.67 


17 


0.13 


0.12 


0.01 


-0.83 


6 


0.09 


0.09 


0.04 


63 077 


-0.27 


2 


0.07 


0.14 


0.04 


-0.67 


11 


0.07 


0.10 


0.02 


-0.82 


5 


0.08 


0.08 


0.03 


76 932 


-0.51 


3 


0.17 


0.13 


0.05 


-0.76 


10 


0.07 


0.11 


0.02 


-1.01 


2 


0.06 


0.08 


0.03 


102 365 


-0.28 


4 


0.08 


0.12 


0.04 


-0.28 


19 


0.06 


0.10 


0.02 


-0.37 


6 


0.02 


0.08 


0.03 


128 620 


+0.29 


6 


0.10 


0.11 


0.04 


+0.26 


20 


0.06 


0.10 


0.02 


+0.29 


6 


0.04 


0.07 


0.02 


131 117 


+0.06 


6 


0.08 


0.10 


0.04 


+0.11 


19 


0.06 


0.10 


0.02 


+0.06 


6 


0.04 


0.08 


0.03 


160691 


+0.25 


4 


0.06 


0.12 


0.04 


+0.28 


19 


0.04 


0.10 


0.02 


+0.32 


6 


0.05 


0.05 


0.01 


196378 











0.04 


-0.38 


14 


0.08 


0.11 


0.02 


-0.44 


3 


0.01 


0.08 


0.03 


199 288 


-0.33 


3 


0.02 


0.13 


0.04 


-0.59 


14 


0.06 


0.11 


0.02 


-0.77 


3 


0.02 


0.08 


0.03 


203 608 











0.05 


-0.69 


9 


0.09 


0.11 


0.01 


-0.70 


3 


0.06 


0.08 


0.03 



Table 9. The same as Table|8] but for Co, Ni, Ce, Nd, and Sm. 



HD 


[Co/H] 


N 


0"[Co/HJ 


uncert.[Co/Hi 


uncert.[Co/Fej 


[Ni/H] 


N 


O" [Ni/H] 


uncert.[Ni/Hi 


uncert.[Ni/Pi.j 


[Ce/H] 


N 


0"[Ce/Hl 


uncert.ice/Hi 


uncert.[c 


2151 


+0.04 


7 


0.03 


0.09 


0.03 


-0.04 


11 


0.05 


0.11 


0.02 


-0.12 


5 


0.10 


0.10 


0.04 


9562 


+0.18 


7 


0.07 


0.09 


0.04 


+0.20 


11 


0.09 


0.10 


0.02 


-0.04 


5 


0.05 


0.10 


0.04 


16417 


+0.16 


8 


0.04 


0.09 


0.04 


+0.14 


11 


0.03 


0.10 


0.02 


+0.03 


5 


0.06 


0.10 


0.04 


20766 


-0.22 


8 


0.06 


0.10 


0.04 


-0.22 


11 


0.07 


0.11 


0.03 


-0.22 


4 


0.06 


0.11 


0.04 


20 807 


-0.20 


8 


0.06 


0.09 


0.04 


-0.23 


10 


0.07 


0.11 


0.03 


-0.31 


4 


0.07 


0.10 


0.04 


22484 


-0.01 


8 


0.08 


0.05 


0.02 


-0.07 


9 


0.05 


0.10 


0.01 


-0.14 


5 


0.07 


0.07 


0.0 


22 879 


-0.55 


6 


0.09 


0.12 


0.04 


-0.79 


9 


0.14 


0.13 


0.04 


-0.72 


5 


0.05 


0.13 


0.05 


30 562 


+0.24 


8 


0.06 


0.09 


0.03 


+0.23 


11 


0.07 


0.10 


0.02 


+0.14 


5 


0.05 


0.10 


0.04 


43 947 


-0.29 


8 


0.06 


0.09 


0.04 


-0.31 


12 


0.10 


0.11 


0.03 


-0.29 


5 


0.07 


0.11 


0.04 


52 298 


-0.23 


5 


0.08 


0.09 


0.03 


-0.38 


5 


0.05 


0.11 


0.03 


-0.27 


4 


0.09 


0.10 


0.04 


59 984 


-0.54 


6 


0.10 


0.13 


0.05 


-0.71 


6 


0.10 


0.14 


0.04 


-0.81 


5 


0.04 


0.14 


0.05 


63 077 


-0.53 


6 


0.07 


0.12 


0.04 


-0.62 


7 


0.09 


0.12 


0.05 


-0.80 


5 


0.03 


0.13 


0.05 


76 932 


-0.61 


4 


0.02 


0.11 


0.04 


-0.80 


7 


0.09 


0.12 


0.04 


-0.81 


5 


0.08 


0.13 


0.05 


102 365 


-0.23 


8 


0.05 


0.10 


0.04 


-0.29 


10 


0.05 


0.11 


0.03 


-0.38 


4 


0.08 


0.11 


0.04 


128 620 


+0.30 


8 


0.06 


0.09 


0.04 


+0.28 


9 


0.04 


0.10 


0.02 


+0.15 


5 


0.08 


0.10 


0.04 


131 117 


+0.05 


8 


0.04 


0.08 


0.03 


+0.12 


12 


0.04 


0.11 


0.02 


-0.03 


5 


0.05 


0.09 


0.04 


160691 


+0.29 


8 


0.04 


0.09 


0.04 


+0.34 


11 


0.04 


0.09 


0.02 


+0.18 


5 


0.06 


0.09 


0.04 


196378 


-0.31 


2 


0.17 


0.09 


0.04 


-0.40 


8 


0.11 


0.12 


0.03 


-0.41 


3 


0.03 


0.11 


0.04 


199 288 


-0.39 


4 


0.05 


0.11 


0.04 


-0.63 


6 


0.05 


0.12 


0.04 


-0.66 


3 


0.07 


0.12 


0.05 


203 608 


-0.18 


3 


0.09 


0.11 


0.04 


-0.74 


5 


0.09 


0.12 


0.03 


-0.73 


4 


0.09 


0.12 


0.05 



HD 


[Nd/H] 


N 


uncert.[Nd/H] 


uncert.[Nd/Fcj 


[Sm/H] 


N uncert.[Sm/Hj 


uncert.[Sm/Fe] 


2151 


-0.03 


2 


0.11 


0.04 


-0.24 


1 0.11 


0.04 


9562 


+0.06 


2 


0.12 


0.04 


+0.06 


1 0.11 


0.04 


16417 


+0.06 


2 


0.12 


0.04 


+0.06 


1 0.11 


0.04 


20 766 


-0.10 


2 


0.12 


0.05 


-0.11 


1 0.12 


0.04 


20 807 


-0.12 


2 


0.12 


0.05 


-0.19 


1 0.11 


0.04 


22484 


+0.04 


2 


0.08 


0.04 


-0.13 


1 0.07 


0.04 


22 879 


-0.29 


2 


0.14 


0.06 


-0.49 


1 0.13 


0.05 


30 562 


+0.16 


2 


0.11 


0.04 


+0.14 


1 0.11 


0.04 


43 947 


-0.17 


2 


0.12 


0.05 


-0.28 


1 0.11 


0.04 


52 298 


-0.10 


2 


0.12 


0.05 


-0.29 


1 0.11 


0.04 


59 984 


-0.47 


2 


0.15 


0.06 


-0.74 


1 0.14 


0.05 


63 077 


-0.40 


2 


0.14 


0.06 


-0.85 


1 0.13 


0.05 


76 932 


-0.44 


2 


0.14 


0.06 


-0.93 


1 0.13 


0.05 


102 365 


-0.13 


2 


0.13 


0.05 


-0.23 


1 0.12 


0.05 


128 620 


+0.28 


2 


0.12 


0.04 


+0.26 


1 0.11 


0.04 


131 117 


+0.10 


2 


0.11 


0.04 


-0.11 


1 0.10 


0.04 


160691 


+0.14 


2 


0.12 


0.03 


+0.28 


1 0.12 


0.04 


196378 











-0.24 


1 0.11 


0.05 


199 288 


-0.32 


1 


0.13 


0.05 


-0.39 


1 0.13 


0.05 


203 608 











-0.62 


1 0.12 


0.05 
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Fig. 11. Percent differences between the EWs of HD 22 879 and 
HD 76 932. Solid lines correspond to mean differences, and the 
dotted ones represent the mean differences modified by their 
standard deviations (±1 cr). 

poor stars, the influence of the continuum placement is stronger 
than that estimated by comparing solar EWs, because EWs are 
considerably lower. On the other hand, it is easier to find good 
continuum windows, because they are wider and more numer- 
ous. So, we compared the EWs of HD 22 879 and HD 76932, 
which are metal-poor stars with very similar atmospheric pa- 
rameters (specially their effective temperatures, which are vir- 
tually identical). The high similarity of atmospheric parameters 
mean that differences in EWs are caused almost exclusively 
by the uncertainties in the EW measurements and continuum 
placement. Only Fe i lines were compared, in order to minimize 
the effect of the stars having slightly different surface gravities. 
The EW percent difference between the two stars present 19% 
dispersion for lines with EW < 15 mA, and 5% for those with 
EW > 15 mA (Fig.ll 1>: the Fe EWs of the metal-poor standard 
stars were modified by these values. 

Log gf are also a source of uncertainty. But, since they are 
derived from solar EWs, their contribution is the same as that 
from the EWs themselves. That is to say, in practice we take 
into account the influence of EWs twice, so as to account for 
the solar EWs per se, and for the log gf values. 

The influence of metallicity in the Fe abundance uncer- 
tainty is the uncertainty itself. This means that the Fe abun- 
dance uncertainty determination is, rigorously speaking, an it- 
erative process. As a first step, Fe abundances were modified 
by the average dispersion between lines of all stars (0.07 dex). 
As can be seen in Table 1 101 the influence of metallicity is very 
weak. The contribution of this source would have to be 0.04 dex 
higher for the total uncertainty to be raised just 0.01 dex. As a 
second step, we should have modified the Fe abundances by the 
total Fe uncertainty just obtained, on the first step, by the RMS 
of all sources (including the metallicity itself). But, given the 
low contribution of this source, we stopped at the first step. 

The total Fe uncertainties for the standard stars were ob- 
tained by the RMS of the six contributions: effective tem- 
peratures, surface gravities, microturbulence velocities, EWs, 
\oggf, and metallicities. Total uncertainties for the other stars 
were obtained by weighted average of the standard stars values, 



using as weight the reciprocal of the distance of the star to each 
standard star, in the T^s vs. [Fe/H] plane. In Fig. ^| we see an 
example, in which the star at the center uses the reciprocal of 
the lengths of the dotted lines as weights. 

Uncertainties for the abundances of the other elements were 
determined the same way as for Fe. Although we determined 
individual [Fe/H] uncertainties for each star, we used a sole 
value of 0.10 dex when estimating the metallicity contribution 
to the uncertainties of the other elements. This can be justi- 
fied by the low dispersion between [Fe/H] uncertainties, and 
by the very low sensitivity of the abundances to the metallicity 
of the model atmosphere. TablellOloresents contributions from 
the individual sources, total uncertainties, and some averages 
for the standard stars, relative to H. Values relative to Fe can be 
obtained simply by subtracting the [Fe/H] value from the [ele- 
ment/H] value (e.g., uncert.[v/Fe] = uncert.[v/H] - uncert.[Fe/H])- 

4.3. Adopted abundances 

Figure [2lpresents the abundance patterns for all contaminat- 
ing elements. Error bars were taken as the dispersion between 
the lines or the estimated [element/Fe] uncertainty, whichever 
is the largest (see Tables|S]and|5). The presence of a few (~ 4%) 
outliers can be noticed. For these stars, these abundances will 
not be used for the Eu and Th spectral syntheses. In these cases, 
the discrepant abundances will be substituted by the value they 
would have if they agreed with an exponential fit to the well 
behaved data points. One might wonder if the discrepant abun- 
dances would be a hint of chemical peculiarity. If this was 
true, one such star would exhibit peculiar abundances for all 
elements produced by one specific nucleosynthesis process, 
and this behaviour is not observed. Mark that the three stars 
which present discrepant abundances (HD 63 077, HD 76932, 
and HD 203 608) are all metal-poor, with [Fe/H] < -0.67, 
and therefore more sensitive to noise and continuum place- 
ment errors. The Nd and V abundances of HD 196 378 and 
HD 203 608, which have no measured lines, were established 
this same way. 

5. Eu and Th abundance determination 

The Eu and Th abundance determination was carried out, for 
our sample stars, using a single absorption line for each el- 
ement. The only acceptably strong Th line, which allows an 
accurate abundance determination, is the Thii line located at 
4019.13 A. This line is, nevertheless, blended with many oth- 
ers, rendering imperative the use of spectral synthesis. One 
other line is often cited in the literature: a Thii line located 
at 4086.52 A. Unfortunately, although less contaminated than 
the line at 4019.13 A, it is too weak to be used. 

Eu has only one adequately strong and uncontaminated 
line, located at 4129.72 A, available for abundance determi- 
nations. However, this fine has a significantly non-Gaussian 
profile by virtue of its HFS and isotope shift. Therefore, spec- 
tral synthesis is also required for this element. In stars of other 
spectral types and other luminosity classes (e.g. giants), it is 
possible to use other Th and Eu lines. This is, nonetheless, not 
true for our sample stars (late-F and G dwarfs/subgiants). 



Table 10. [element/H] abundance uncertainties for die standard stars. Last column contains the average uncertainty for each row, and the last row contains 
the average total uncertainty for the four standard stars. 



Parameter 


AParameter 


HD 


[Fe/H] 


[Ti/H] 


[V/H] 


[Cr/H] 


[Mn/H] 


[Co/H] 


[Ni/H] 


[Ce/H] 


[Nd/H] 


[Sm/H] 


Average 






160691 


+0.01 


+0.02 


+0.03 


+0.02 


+0.02 


+0.02 


+0.02 


+0.00 


+0.01 


+0.01 


+0.02 ±0.01 


T a 


+27 K 


22484 


+0.02 


+0.02 


+0.03 


+0.02 


+0.03 


+0.03 


+0.02 


+0.01 


+0.01 


+0.01 


+0.02 ±0.01 


^ eft 


63 077 


+0.02 


+0.02 


+0.03 


+0.01 


+0.02 


+0.02 


+0.02 


+0.01 


+0.01 


+0.01 


+0.02 ±0.01 






59 984 


+0.02 


+0.02 


+0.03 


+0.02 


+0.03 


+0.02 


+0.02 


+0.01 


+0.02 


+0.01 


+0.02 ±0.01 






160691 


+0.00 


+0.00 


+0.00 


+0.00 


+0.00 


+0.00 


+0.00 


+0.01 


+0.01 


+0.01 


+0.00 ±0.00 




+0.02 dex 


22484 


+0.00 


+0.00 


+0.01 


+0.00 


+0.00 


+0.01 


+0.00 


+0.01 


+0.01 


+0.01 


+0.01 ±0.01 


63 077 


+0.00 


+0.00 


+0.00 


+0.01 


+0.00 


+0.00 


+0.00 


+0.01 


+0.01 


+0.01 


+0.00 ±0.01 






59 984 


+0.00 


+0.01 


+0.00 


+0.00 


+0.00 


+0.02 


+0.00 


+0.01 


+0.01 


+0.01 


+0.01 ±0.01 




+0.05 km 


160691 


-0.01 


-0.02 


+0.00 


-0.01 


-0.01 


+0.00 


-0.02 


-0.01 


-0.01 


-0.01 


-0.01 ±0.01 


? 


22484 


-0.01 


-0.01 


-0.01 


-0.02 


+0.00 


-0.01 


-0.01 


-0.01 


-0.01 


-0.01 


-0.01 ±0.00 


+0.23 km s"' 


63 077 


-0.03 


-0.03 


-0.01 


-0.04 


+0.00 


+0.00 


+0.00 


-0.02 


-0.01 


-0.01 


-0.01 ±0.01 




59 984 


-0.03 


-0.03 


+0.00 


-0.04 


+0.00 


+0.00 


-0.01 


-0.01 


-0.01 


-0.01 


-0.01 ±0.01 


EW 
and 
logg/ 


EW<20mA:+15% 


160691 


+0.06 


+0.07 


+0.04 


+0.07 


+0.06 


+0.03 


+0.07 


+0.04 


+0.05 


+0.04 


+0.05 ±0.01 


EW > 20 mA: +5% 


22484 


+0.06 


+0.06 


+0.08 


+0.07 


+0.03 


+0.06 


+0.06 


+0.06 


+0.08 


+0.08 


+0.06 ±0.02 


EW < 15 mA: +19% 


63 077 


+0.07 


+0.05 


+0.09 


+0.08 


+0.06 


+0.09 


+0.10 


+0.10 


+0.10 


+0.10 


+0.08 ±0.02 


EW > 15 mA: +5% 


59 984 


+0.06 


+0.06 


+0.10 


+0.06 


+0.05 


+0.08 


+0.08 


+0.09 


+0.10 


+0.09 


+0.08 ±0.02 


[Fe/H] 


Fe: +0.07 dex 


160691 

/^i^ AO A 

22484 


+0.01 
+0.00 


+0.01 
+0.01 


+0.01 
+0.00 


+0.01 
+0.01 


+0.00 
+0.00 


+0.01 
+0.01 


+0.01 
+0.01 


+0.03 
+0.03 


+0.04 
+0.03 


+0.04 
+0.03 


+0.02 ±0.01 
+0.01 ±0.01 


OtVlPf plplTlPTltQ' 


OiUl / 


+0.0U 


+0.01 


+0.01 


+0.U1 


1 A An 
+0.0U 


1 A AA 
+ 0.1)0 


1 A AA 


1 A AT 
+U.UI 


1 A AQ 
+ 0.03 


1 A A 1 
+U.U1 


+0.01 ±0.01 




_i_n 1 n H<=»v 


59 984 


+0.00 


+0.01 


+0.00 


+0.01 


+0.01 


+0.01 


+0.01 


+0.02 


+0.02 


+0.01 


+0.01 ±0.01 






160691 


0.09 


0.09 


0.12 


0.10 


0.05 


0.09 


0.09 


0.09 


0.12 


0.12 


0.08 ±0.02 


Total 




22484 


0.09 


0.10 


0.06 


0.10 


0.09 


0.05 


0.10 


0.07 


0.08 


0.07 


0.10 ±0.02 


uncertainty 




63 077 


0.09 


0.09 


0.14 


0.10 


0.08 


0.12 


0.12 


0.13 


0.14 


0.13 


0.12 ±0.02 






59984 


0.11 


0.08 


0.13 


0.12 


0.09 


0.13 


0.14 


0.14 


0.15 


0.14 


0.11 ±0.02 


Average 






0.10 


0.09 


0.11 


0.11 


0.08 


0.10 


0.11 


0.11 


0.12 


0.12 


0.10 ±0.03 


total 






±0.01 


±0.01 


±0.04 


±0.01 


±0.02 


±0.04 


±0.02 


±0.03 


±0.03 


±0.03 



Note: Contributions from EW and log gf are equal. The values which are presented correspond to each contribution, and not to the total of these two sources. 
Therefore, they are taken in consideration twice when calculating the total uncertainties. 
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The code used to calculate the synthetic spectra was kindly 
made available to us by its developer, Monique Spite. This code 
calculates a synthetic spectrum based on a list of lines with their 
respective atomic parameters (log gf, and excitation potential 
of the lower level of the transition in eV). We employed the 
same model atmospheres and partition functions used in the de- 
termination of atmospheric parameters and abundances of con- 
taminating elements (Section f3.U . The only exception was for 
Th, whose partition functions were calculated by us, based on 
data from H. Holweger that was published as a private commu- 
nication in MKB92. Irwin's data are too inaccurate, because 
for the singly and doubly ionised states of Th, he scales the 
partition functions along iso-electronic sequences from a lower 
mass element, due to lack of specific Th data. This resulted in 
values much lower than those of Holweger, which were calcu- 
lated by summing a large number of Th atomic levels. 

5.1. Europium 

5.1 .1 . Spectral synthesis 

The list of lines used for calculation of the synthetic spectrum 
of the Eu region was based on the line lists of th e Vienna 
Atomic Line Database - VALD (Kuoka et al. '1999^. Initially, 
we took all Hnes found between 4128.4 A and 4130.4 A. The 
first cut eliminated all lines of elements ionised two or more 
times, and all lines with lower-level excitation potential greater 
than 10 eV, since these lines are unobservable in the photo- 
spheres of stars with effective temperatures as low as those 
in our sample. As a next step, we calculated synthetic spec- 
tra for the four stars with the most extreme effective tempera- 
tures and metallicities in our sample (HD 22484, HD 59984, 
HD 63 077, and HD 203 608) using laboratory log gf values, 
and the atmospheric parameters and abundances previously 
obtained. Lines with EWs lower than 0.01 mA in these four 
stars were removed. T he Sc i line located at 4129.750 A, cited 
as important by Mashonkina & Gehren ( 2000 ), presents neg- 
ligible EW in all four standard stars. 23 Ti, V, Cr, Fe, Co, 
Nb, Ce, Pr, Nd, Sm, and Dy lines were kept, besides the Eu 
line itself. Adopted wavelengths were taken from the VALp 
list llRardfcKocklll994 lKunic7lll993[ 1 1994 IWhalindfl983t 
IWickUffe et alJll994 . To improve the fits further, three artifi- 
cial Fei lines were added as substitutes for unkn own blends 
jLawleret al.ll200lt IMashonkina & GehrenlEoOOl) : the influ- 
ence of these artificial lines in the obtained Eu abundances is 
nevertheless small, as they are located relatively far from the 
center of the Euii line. The final list is presented in Tablell II 

In the spectral synthesis calculations, the Euii line 
was substituted by its HFS components, ca lculated by 
us using Casimir's e quation dCasimid Il963l). and HPS 
labora t ory data f rom Becker etal. (1993), Brostrom et al. 
( 1 19951) . IViUemoes & Wang (199 4), and MoUer et al. (1993). 
Isotope shift was t aken into account, using data from 
iBro strom et al.l ( 1 19951) and the solar abundance isotopic ratio 
e(i5iEu)/e('53Eu) = 1.00 ±0.29 O^awler et al. 2001 ). The solai- 
value was used for all stars because there are strong indications 
that it is not metallicity-dependent, since several works dealing 
with r-process nucleosynthesis in very-metal-poor halo stars 



Cr I and Cell Sm II Nb I 




I I I I I I I I I I I I I I I : : I 

4128.8 4129.2 4129.6 4130.0 

X(k) 

Fig. 13. Synthesized spectral region for the Eu ii line at 
4129.72 A, for the Sun. The two thick lines represent the total 
synthetic spectrum and the Eun line. The thin lines represent 
the most dominant lines in the region, including two artificial 
Fe I lines. Points are the observed FEROS solar spectrum. 

find simila r results ('Pfeiffer et al.' 'l997l ICowan et alJ 1 19991: 

ISneden et al 2002: Aoki et al...2003.) . 

We adopted a projected rotation ve locity for the inte- 
grated solar disk v sin ; = 1.8 km s"' (Mashonkina '2000'; 
IMashonkina & GehrenI l2000l) . The width of the Gaussian 
profile convolved to the synthetic spectrum to take macro- 
turbulence and instrumental broadenings into account was 
obtained by fitting the Fen located at 4 128.742 A 
which is sufficiently isolated for this purpose (IWoolfetalJ 
199?, iKoch & Edva rdsson 2002). Solar log gf were ob- 
tained by keeping abun dances fixed at their solar values 
dGrevesse & Sauval 19981) . and fitting the observed solar spec- 
trum. Only the stronger Ti, Cr, and Fe lines had their log g f 
adjusted. For the Eu n line, we adopted a fixed log gf from 
iKomarovskiil (Il99ll) : this value was distributed among the 
HFS components according to relative intensities calculated 
by White & Eliason, tabulated in Condon & Shortlev C 19671) . 
The other, weaker lines had their logef kept at their labora- 
tory v alues adopted from the VALD fist ( Kurucz 19931 1 1994 
IWickliffe et aL. 1994 ). As we kept the Eu log gf fixed, we al- 
lowed its abundance to vary. The complete solar spectrum fit 
procedure was carried out iteratively because it requires al- 
ternate Gaussian profile adjustments (which modify the shape 
of the lines) and log gf adjustments (which modify the EWs 
of the lines). The solar log gf and Gaussian profile determi- 
nation was accomplished independently for the two sources 
of Eu spectra used, leading to different Eu solar abundances 
(loge(Eu)o = -hO.37 and -hO.47 for the CES and FEROS, re- 
spectively). The stellar [Eu/H] abundance ratios, used in this 
work in lieu of the absolute abundances, were obtained by sub- 
tracting the appropriate solar value from the stellar absolute 
abundances ([Eu/H] - loge(Eu)star - loge(Eu)0). Figure [Oj 
presents the observed and synthetic solar spectrum in the Eu re- 
gion. The strongest lines that compose the total synthetic spec- 
trum are presented independently. 
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Table 11. Line list used in the spectral synthesis of the Eu ii line at 4129.72 A. 
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A (A) 


Element 


r (eV) 


loe S frp^ 


l02 S fpPTin^ 


\o^gf source 


Obs. 


4128.742 


Fe II 


2.580 


-3.832 


-3.554 


solar tit 




4129.000 


Cri 


4.210 


-2.603 


-2.603 


VALD 




4129.040 


Rui 


1.730 


-1.030 


-1.030 


VALD 




4129.147 


Prii 


1.040 


-0.100 


-0.100 


VALD 




4129.159 


Tin 


1.890 


-2.330 


-2.210 


solar fit 




4129.166 


Til 


2.320 


+0.131 


+0.251 


solar fit 




4129.165 


Cri 


3.010 


-1.948 


-1.948 


VALD 




4129.174 


Ceil 


0.740 


-0.901 


-0.901 


VALD 




4129.196 


Cri 


2.910 


-1.374 


-1.254 


solar fit 




4129.220 


Fe I 


3.420 


-2.280 


-2.160 


solar fit 




4129.220 


Smii 


0.250 


-1.123 


-1.123 


VALD 




4129.425 


Dy II 


0.540 


-0.522 


-0.522 


VALD 




4129.426 


Nbi 


0.090 


-0.780 


-0.780 


VALD 




4129.461 


Fe I 


3.400 


-2.180 


-1.920 


solar fit 




4129.530 


Fe I 


3.140 


-3.425 


-3.455 


solar fit 


artificial 


4129.610 


Fe I 


3.500 


-3.700 


-3.730 


solar fit 


artificial 


4129.643 


Til 


2.240 


-1.987 


-1.987 


VALD 




4129.657 


Til 


2.780 


-2.297 


-2.297 


VALD 




4129.721 


Eu II 


0.000 


+0.173 


+0.173 


Komarovskii (199D 


HFS 


4129.817 


Co I 


3.810 


-1.808 


-1.808 


VALD 




4129.837 


Ndii 


2.020 


-0.543 


-0.543 


VALD 




4129.965 


Fei 


2.670 


-3.390 


-3.290 


solar fit 


artificial 


4129.994 


Vi 


2.260 


-1.769 


-1.769 


VALD 




4130.037 


Fei 


1.560 


-4.195 


-4.030 


solar fit 




4130.038 


Fei 


3.110 


-2.385 


-2.225 


solar fit 




4130.068 


Fei 


3.700 


-3.763 


-3.763 


VALD 




4130.073 


Cri 


2.910 


-1.971 


-1.971 


VALD 





References: see text. 



For the Eu synthesis in the sample stars, we kept the abun- 
dances of all elements other than Eu fixed at the values deter- 
mined using EWs, in Sect.|3] Small adjustments were allowed 
only to the abundances of Ti, Cr, and Fe, in order to improve 
the fit. The adjustments were kept within the uncertainties of 
the abundances of these elements - see Sect. |3 The abun- 
dances of elements for which we did not measure any EW were 
determined scaling the solar abundances, following elements 
produced by the same nucleosynthetic processes. We chose to 
scale the Sm abundances as well, even though they were deter- 
mined by us, because our results for this element exhibit high 
scatter and absence of well-defined behaviour (see Fig. I12> . 
Sm and Dy abundances were obta ined followin g Eu (produced 
mainly by r-process, according to lBurris et al.l t2000): Nb, Rb, 
and Pr were obtained following Nd (produced by both r-process 
and s-process, also according to Buii'is et al. 2000) . Mark that 
the influence of these lines is marginal, since they are weak 
and are not close to the Eu line. Projected rotation veloci- 
ties V sin i and macroturbulence and instrumental broadening 
Gaussian profiles were determined by fitting the Fe ii line at 
4128.742 A, as done with the solar spectrum. 

Stars observed more than once had their spectra analysed 
independently, and the results were averaged. Abundances ob- 
tained with the different spectra of one object presented a maxi- 
mum variation of 0.02 dex. Examples of spectral syntheses can 
be seen in Fig.[2]for the two stars with extreme Eu abundances 
that were observed with both FEROS and CES (HD 63 077 and 
HD 160691). Note the effect of CES spectra higher sampling. 



which better puts the asymmetric profile of the Eu ii line in ev- 
idence. 

Eu abundances obtained using CES and FEROS spectra are 
compared in Fig. ^] A linear fit was calculated, and the re- 
sulting relation was used to convert FEROS results into the 
CES system. This way, we ended up with two sets of Eu abun- 
dances: one totally homogeneous, obtained with FEROS spec- 
tra, and one obtained partially with CES spectra and partially 
with corrected-into-CES-system FEROS spectra. We plotted a 
[Eu/Fe] vs. [Fe/H] diagram for both sets of abundances, pre- 
sented in Fig.^] Linear fits calculated for each set show that 
the CES abundances present a lower scatter cr, and for this rea- 
son the CES results were adopted for all subsequent analyses. 

5.1.2. Uncertainty assessment 

[Eu/H] abundance uncertainties were estimated with a proce- 
dure identical to the one used for the abundances of contam- 
inating elements determined with EWs, thoroughly described 
in Sect. 14.21 It consists basically of varying the atmospheric pa- 
rameters and continuum position independently, by an amount 
equal to their uncertainties, and recalculating the synthetic 
spectra. Total uncertainties are obtained by RMS of the individ- 
ual sources of uncertainty. We chose four stars as standards for 
the assessment, with extreme eff'ective temperatures and metal- 
licities: HD 160691 (cool and metal-rich), HD 22484 (hot and 
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Fig. 15. Comparison between Eu abundances determined using 
CES and FEROS spectra. Solid line is the identity. Dotted line 
is a linear fit, whose parameters are shown in the figure. 
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Fig. 16. [Eu/Fe] vs. [Fe/H] diagram for all sample stars. 
Abundances obtained with CES and FEROS spectra are pre- 
sented. The Sun is not included in the comparison, since its 
abundance is always equal to zero by definition. Solid and dot- 
ted lines are linear fits to these two data sets, whose scatters are 
shown in the figure. Mark that data from the CES are slightly 
better than those from FEROS. 



continuum placement clearly dominate. Microturbulence ve- 
locity variations have no effect in the [Eu/H] abundances. This 
happens because the HFS broadens the l ine so much that it 



122^ 



mm 



becomes nearl y unsatu rated, as found by IWoolf et a^J 
WTL95), Lawl eret all dioOl) and Koch & Edvardsso^ 
KE02). [Element/Fe] uncertainties were obtained simply by 
subtracting the [Fe/H] uncertainty from the [element/H] value 
(e.g., uncert. [Eu/Fe] = uncert.[Eu/H] - uncert.^e/H])- Total uncer- 
tainties for the other sample stars were obtained by weighted 
average of the standard stars values, using as weight the recip- 
rocal of the distance of the star to each standard star, in the Tgff 
vs. [Fe/H] plane; this procedure is identical to the one used 
for the uncertainties of contaminating elements (Sect. 14. 2> . 
Table [Hpresents the final [Eu/H], [Th/H], and [Th/Eu] abun- 
dance ratios, for all sample stars, along with their respective 
uncertainties relative to H and Fe. 

Non-LTE effects in the abundances obtained with 
the Euii line at 4129.72 A have been calculated by 
'Mashonkina & Gehren (2000). In the line formation layers, 
the ground state was found to be slightly underpopulated, and 
the excited level overpopulated. As a consequence, our results 
would show a small difference if our method employed abso- 
lute abundances. But since we employ a differential analysis, 
these differences cancel each other partially, becoming negligi- 
ble (WTL95 and KE02). 



5.1.3. Comparison with literature results 



We conducted a painstaking search of the literature for works 
with Eu abundances, and selected WTL95 and KE02 as the 
ones with the most careful determinations of disk stars, as 
well as a sizable sample of high statistical significance . Thei r 
sample is composed of a subset from Edvardsson et a l1 (ll99:^ . 
and they use the same atmospheric parameters (obtained from 
Stromgren photometric calibrations). The Eu abundances were 
determined by spectral synthesis using the same line we used, 
with a procedure fundamentally identical to ours. KE02 merged 
their database with that from WTL95 by means of a simple 
linear conversion, obtained by intercomparison. We used this 
merged set of abundances to compare to our results. 



metal-rich), HD 63 077 (cool and metal-poor), and HD 203 608 
(hot and metal-poor) - Fig.FTOI 

Continuum placement is the preponderant source of Eu 
abundance uncertainty. The placement uncertainty itself is 
higher for metal-rich stars, as these present stronger absorp- 
tion lines, lowering the apparent continuum which restricts the 
choice of good normalisation windows. To estimate the influ- 
ence of continuum placement uncertainties, we multiplied the 
spectral flux of the metal-rich standard stars by 0.98, and that 
of the metal-poor standard stars by 0.99, and recalculated the 
synthetic spectra. 

Tableinipresents the [Eu/H], [Eu/Fe], [Th/H], [Th/Fe], and 
[Th/Eu] uncertainties for each standard star. Metallicity and 



FigureFTTlpresents a [Eu/H] vs. [Fe/H] diagram with our re- 
sults and those from WTL95 and KE02. Both data sets exhibit 
the same behaviour, but our results present scatter 36% lower 
than those of WTL95 and KE02, even though our sample is 
6 times smaller than theirs (linear fits result in 0-0^. ^g^]^ - 0.050 
and crwTL95/KE02 - 0.078). The lower scatter of our results 
is a consequence of improvements introduced in our analysis, 
like the use of atmospheric parameters obtained by us through 
the detailed and totally self-consistent spectroscopic analysis 
carried out in Sect. |3] WTL95 and KE02 demonstrated, using 
Monte Carlo simulations of observational errors, that the ob- 
served scatter is not real, but mainly a result of observational, 
analytical, and systematic errors. 
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Fig. 17. [Eu/H] vs. [Fe/H] diagram for our sample stars and 
those from WTL95 and KE02. Average error bars for the two 
data sets are provided in the lower right corner. Note that the 
behaviour of both data sets is similar, but that our abundances 
present considerably lower scatter. 



5.2. Thorium 

5.2.1. Spectral synthesis 

The list of lines used for calculation of the synthetic spectrum 
of the Th region was constructed following the same procedure 
used for the Eu analysis. 14 V, Mn, Fe, Co, Ni, W, Ce, and U 
lines were kept, besides the Th line itself. Adopted wavelengths 
were taken from the VALD list (Bard & Koc k .1994; Kuruc?^ 
Il993lll994trwhaling|l983'l) . with the exception of the Fe i Une 
at 4019.043 A, Nil at 401 9.067 A, a nd Thii at 4019.130 A, 
determined in laboratory bv lLearneret al . (1991) with a preci- 
sion higher than 0.0005 A. The Co i lines were substituted by 
their HFSs components, calculated by us using Casimir's equa- 
tion, and HFS laboratory dat a from Chil d s & Goodrnan ( 1 968). 
Guthohrlein & KelleJ (Il990l) . IPickering & Thornd ri99( 



and 



Pickering ' ('199^ To improve the fits further, one artificial Fe ii 



line was added as substitute for unknown blends (MKB92; 
iFranqois et aLlil993l) : the influence of this artificial line in the 
obtained Th abundances is nevertheless small, as it is located 
relatively far from the center of the Thii line. The final list is 
presented in Tablell2l 

The width of the Gaussian profile used to take macroturbu- 
lence and instrumental broadenings into account was obtained 
by fitting the Co i fines located at 4019.293 A and 4019.297 A. 
Solar log gf were obt ained by keeping abundances fixed at 
their solar values ( Grevesse & Sauval 199 8). and fitting the ob- 
served solar spectrum. Only the stronger Mn, Fe, Ce, and Ni 
lines had their log gf adjusted. For the Thii line, we adopted 
a fixed log gf from J^Iilssqn et al... C2002.): for the Co i lines we 
adopted fixed log g f from lLawler et al.ll99(il r4019 1 32 A) and 
VALD (4019.293 A and 4019.297 A, 'Kurucz'l994'). The other, 
weaker lines had their log g f kept at their labor atory values 
adopted from flie VALD list teuruczll993lll994 . As we kept 
the Th log gf fixed, we allowed its abundance to vary. The solar 
log gf and Gaussian profile determination was accompfished 
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Fig. 18. Synthesized spectral region for the Th ii line at 
4019.13 A, for the Sun. The two thick lines represent the total 
synthetic spectrum and the Thii line. The thin lines represent 
the most dominant lines in the region, including an artificial 
Fe II line. Points are the observed CES solar spectrum, obtained 
with the 3.60 m telescope. 



independently for each telescope used, leading to slightly dif- 
ferent Th solar abundances (loge(Th)0 = -1-0.03 and -h0.04 
for the 3.60 m and CAT, respectively). As with Eu, the stel- 
lar [Th/H] abundance ratios were obtained by subtracting the 
appropriate solar value from the stellar absolute abundances 
([Th/H] = loge(Th)star - loge(Th)0). Figure^] presents the 
observed and synthetic solar spectrum in the Th region. The 
strongest lines that compose the total synthetic spectrum are 
presented independently. 

For the Th synthesis in the sample stars, we kept the abun- 
dances of all elements other than Th fixed at the values deter- 
mined using EWs, in Sect.|3 Small adjustments were allowed 
only to the abundances of Mn, Fe, Ce, and Ni, in order to im- 
prove the fit. The adjustments were kept within the uncertain- 
ties of the abundances of these elements - see Sect. ^ The W 
abundances, for which we did not measure any EW, were de- 
termined scaling the solar abundances, following the Fe abun- 
dances. The U abundances were estimated from a simplifying 
hypothesis that the [U/H] stellar abundance at the time of for- 
mation is independent of metallicity. We allied this hypothesis 
to a simple linear age-metallicity relation, in which a star with 
[Fe/H] = -1.00 is 10 Gyr old, and one with [Fe/H] = H-0.00 
is 4.57 Gyr old, i.e., age(Gyr) = 4.57 - 5.43 [Fe/H]. Given 
the known ^^^U half-life t\/2 - 4.46 Gyr, we reach a relation 
[U/H] = -0.50 4- 0.37 [Fe/H]. Mind that the W and U fines 
have only a small influence on the final Th derived by the anal- 
ysis, since they are weak (with EW < 0.2 mA for the metal- 
richest stars) and are not close to the Th line. We employed 
the projected rotation velocities v sin / determined in the Eu 
analysis, and the macroturbulence and instrumental broaden- 
ing Gaussian profiles were obtained by fitting the Co i lines at 
4019.293 A and 4019.297 A, like done wifli die solar spectrum. 

Stars observed more than once had their spectra analysed 
individually, and the results were averaged. Different spec- 



20 E.F. del Peloso et al.: The age of the Galactic thin disk from Th/Eu nucleocosmochronology. I 

Table 12. Line list used in the spectral synthesis of the Thii line at 4019.13 A. 
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Fe II 


3.00 
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solar fit 
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4019.228 


Wi 
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-2.200 


VALD 




4019.293 


Co I 
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-3.232 


-3.232 


VALD 


HPS 


4019.297 


Co I 


0.63 


-3.769 


-3.769 


VALD 


HPS 



References: see text. 



tra for each star presented a 0.02 dex maximum variation. 
Examples of spectral syntheses can be seen in Fig.[^for the 
two stars with extreme Th abundances that were observed with 
both 3.60 m and CAT (HD 20 766 and HD 128 620). Mark that, 
due to problems that degraded the resolving power of spectra 
obtained with the 3.60 m telescope, these are almost identical 
to the ones obtained with the CAT 

Th abundances determined using spectra obtained with the 
3.60 m and CAT are compared in Fig.|20l A linear fit was cal- 
culated, and the resulting relation was used to convert 3.60 m 
results into the CAT system, and vice versa. This way, we ended 
up with two sets of Th abundances: one obtained partially with 
3.60 m spectra and partially with corrected-into-3.60-m-system 
CAT spectra, and one that is just the opposite, obtained par- 
tially with CAT spectra and partially with corrected-into-CAT- 
system 3.60 m spectra. We plotted a [Th/Fe] vs. [Fe/H] diagram 
for both sets of abundances, presented in Fig. |^ Linear fits 
calculated for each set show that the CAT abundances present 
a lower scatter cr, and for this reason the CAT results were 
adopted for all subsequent analyses. Thus, our adopted Eu and 
Th were obtained in one and the same system, related to one 
single instrument (CES) and telescope (CAT), reinforcing the 
homogeneity of the analysis. 
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Fig. 20. Comparison between Th abundances determined using 
CES spectra obtained with the CAT and with the 3.60 m. Solid 
line is the identity. Dotted line is a linear fit, whose parameters 
are shown in the figure. 

that, depending on the metallicity of the star, has an EW com- 
parable to Th itself. 



5.2.2. Uncertainty assessment 

[Th/H] abundance uncertainties were obtained with a proce- 
dure identical to the one used for Eu. Tablell3lcontains the re- 
sults of the Th uncertainty analysis for the four standard stars, 
and Tablell4lcontains the abundances for all sample stars, with 
their respective uncertainties. Dependence of Th abundances 
on atmospheric parameters is similar to Eu, but exhibit 3 to 
5 times higher sensitivity to continuum placement variations. 
This results from the much lower EW of the Th line, which 
ranges from 1 mA to 8 mA for our sample stars, whereas Eu 
ranges from 20 mA to 90 mA. Furthermore, the Eu line is prac- 
tically isolated from contaminations, while the Th line is lo- 
cated in the red wing of a feature composed of strong Mn, Fe, 
Ce, and Ni lines, and is blended with one V i and one Co i line 



5.2.3. Comparison with literature results 

There are only two works available in the literature in which 
Th abundance s have been determined for Galactic disk stars: 
'da Silva et al.' ( Il990l) and MKB92. In the classic work of 
Butcher ( 1987), which was the first to propose the use of stellar 
Th abundances as chronometers, the ratio between the EWs of 
Th and Nd is adopted as an estimate of [Th/Nd], but no inde- 
pendent Th abundanc es are derived . We did not compare our 
results to those from Ida Silva et alJ because only preliminary 
data is presented for a small sample of four stars. 

MKB92 is are-analysis of the .Butcher ( 1987) data, with the 
same stellar sample and spectra, but making use of more up- 
to-date model atmospheres and more detailed spectral synthe- 
ses. Unfortunately, MKB92 do not present a detailed Th abun- 
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Table 13. [Eu/H], [Eu/Fe], [Th/H], [Th/Eu], and [Th/Eu] abundance uncertainties for the standard stars. 
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Parameter 


AParameter 


HD 


[Eu/Hl 


[Eu/Fe] 


TTh/Hl 


iTh/Fel 

L ± 11/1 tj 


[Th/Eu] 






160 691 


+0.01 


+0.00 


+0.00 


-0.01 


-0.01 




+27 K 


22484 


+0.01 


-0.01 


+0.01 


-0.01 


+0.00 




63 077 


+0.01 


-0.01 


+0.01 


-0.01 


+0.00 






203 608 


+0.02 


+0.00 


+0.01 


-0.01 


-0.01 






160 691 


+0.01 


+0.01 


+0.01 


+0.01 


+0.00 




+0.02 dex 


22484 


+0.01 


+0.01 


+0.01 


+0.01 


+0.00 


63 077 


+0.01 


+0.01 


+0.01 


+0.01 


+0.00 






203 608 


+0.01 


+0.01 


+0.01 


+0.01 


+0.00 




+0.05 km 


160 691 


+0.00 


+0.01 


+0.00 


+0.01 


+0.00 




22484 


+0.00 


+0.01 


+0.00 


+0.01 


+0.00 


+0.23 km 


63 077 


+0.00 


+0.03 


+0.00 


+0.03 


+0.00 




203 608 


+0.00 


+0.03 


+0.00 


+0.03 


+0.00 






160 691 


+0.03 


+0.02 


+0.04 


+0.03 


+0.01 


[Fe/H] 


+0.10 dex 


22484 


+0.02 


+0.02 


+0.03 


+0.03 


+0.01 


63 077 


+0.02 


+0.02 


+0.02 


+0.02 


+0.00 






203 608 


+0.01 


+0.01 


+0.02 


+0.02 


+0.01 




+2% 


160 691 


+0.02 


-0.04 


+0.09 


+0.03 


+0.07 


Continuum 


22484 


+0.03 


-0.03 


+0.10 


+0.04 


+0.07 


+ 1% 


63 077 


+0.02 


-0.05 


+0.11 


+0.04 


+0.09 




203 608 


+0.03 


-0.03 


+0.14 


+0.08 


+0.11 






160 691 


0.04 


0.05 


0.10 


0.05 


0.07 


Total 




22484 


0.04 


0.04 


0.11 


0.05 


0.07 


uncertainty 




63 077 


0.03 


0.06 


0.11 


0.06 


0.09 






203 608 


0.04 


0.04 


0.14 


0.09 


0.11 



Table 14. [Eu/H], [Th/H], and [Th/Eu] abundance ratios for all sample stars, with their respective uncertainties relative to H and 
Fe. 



HD 


[Eu/H] 


uncert.[Ei,/Hj 


uncert.[E„/Fc] 


[Th/H] 


uncert.[Th/Hj 


uncert.[Th/Fc] 


[Th/Eu] 


uncert. [Th/Eu] 


2151 


-0.04 


0.04 


0.05 


+0.00 


0.11 


0.06 


+0.04 


0.08 


9562 


+0.05 


0.04 


0.05 


+0.04 


0.11 


0.06 


-0.01 


0.08 


16417 


+0.12 


0.04 


0.05 


+0.06 


0.11 


0.06 


-0.06 


0.08 


20766 


-0.15 


0.04 


0.05 


-0.25 


0.11 


0.06 


-0.10 


0.08 


20 807 


-0.12 


0.04 


0.05 


-0.17 


0.11 


0.06 


-0.05 


0.08 


22484 


-0.01 


0.04 


0.04 


+0.04 


0.11 


0.05 


+0.05 


0.07 


22 879 


-0.42 


0.04 


0.05 


-0.58 


0.12 


0.07 


-0.16 


0.09 


30562 


+0.23 


0.04 


0.05 


+0.19 


0.11 


0.06 


-0.04 


0.08 


43 947 


-0.18 


0.04 


0.05 


-0.11 


0.12 


0.06 


+0.07 


0.08 


52 298 


-0.13 


0.04 


0.04 


-0.09 


0.12 


0.06 


+0.04 


0.09 


59 984 


-0.48 


0.04 


0.05 


-0.42 


0.12 


0.07 


+0.06 


0.09 


63 077 


-0.48 


0.03 


0.06 


-0.48 


0.11 


0.06 


+0.00 


0.09 


76932 


-0.48 


0.04 


0.05 


-0.54 


0.12 


0.07 


-0.06 


0.09 


102 365 


-0.11 


0.04 


0.05 


-0.23 


0.11 


0.06 


-0.12 


0.08 


128 620 


+0.16 


0.04 


0.05 


+0.09 


0.11 


0.06 


-0.07 


0.08 


131 117 


+0.07 


0.04 


0.05 


+0.10 


0.11 


0.06 


+0.03 


0.08 


160 691 


+0.16 


0.04 


0.05 


+0.09 


0.10 


0.05 


-0.06 


0.07 


196 378 


-0.22 


0.04 


0.04 


-0.31 


0.12 


0.07 


-0.10 


0.09 


199 288 


-0.30 


0.03 


0.05 


-0.32 


0.11 


0.06 


-0.02 


0.09 


203 608 


-0.50 


0.04 


0.04 


-0.46 


0.14 


0.09 


+0.04 


0.11 



dance uncertainty analysis. The authors investigated the influ- 
ence of effective temperature, surface gravity and continuum 
placement variations, but overlooked metallicity and microtur- 
bulence velocity. No estimate of total uncertainty is presented. 
Therefore, we decided not to include error bars in Fig. |22] 
which depicts [Th/H] vs. [Fe/H] for our results and those from 
MKB92. The two data sets exhibit a similar behaviour, but our 



results present scatter 61% lower than those of MKB92 (from 
linear fits we get cToui-woi-k = 0.064 and crMKB92 = 0.165). 
The lower scatter of our results is a consequence of many en- 
hancements introduced in our analysis, like the use of homo- 
geneous and precise atmospheric parameters obtained by us, 
while MKB92 used multiple literature sources, and refinements 
in the spectral synthesis. 
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Fig. 21. [Th/Fe] vs. [Fe/H] diagram for all sample stars. 
Abundances determined using CES spectra obtained with the 
CAT and the 3.60 m are presented. The Sun is not included 
in the comparison, since its abundance is always equal to zero 
by definition. Solid and dotted lines are linear fits to these two 
data sets, whose scatters are shown in the figure. Mark that data 
from the CAT are better than those from 3.60 m. 
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Fig. 22. [Th/H] vs. [Fe/H] diagram for our sample stars and 
those from MKB92. Note that the behaviour of both data sets 
is similar, but that our abundances present considerably lower 
scatter. 

6. Conclusions 

Homogenous, fully self-consistent atmospheric parameters, 
and their respective uncertainties, have been obtained for all 
sample stars. The use of two different and internally homoge- 
neous criteria for effective temperature determination allowed 
us to achieve a very low uncertainty (27 K) for this impor- 
tant parameter. Uncertainties of the other parameters were also 
found to be adequate for our needs. 

Abundances of the elements that contaminate the Th n and 
Eu II spectral regions - Ti, V, Cr, Mn, Co, Ni, Ce, Nd and Sm - 
have been determined by detailed spectroscopic analysis, rela- 
tive to the Sun, using EWs. For the elements with meaningful 
HFSs (V, Mn, and Co), this has been taken into account. A 



thorough estimation of the uncertainties has been carried out. 
The average [element/H] uncertainty - (0.10 ± 0.02) dex - was 
found to be satisfactorily low. 

Eu and Th abundances have been determined for all sample 
stars using spectral synthesis of the Euii line at 4129.72 A and 
the Thii line at 4019.13 A. Comparison of our results with the 
literature shows that our analysis yielded a similar behaviour, 
but with considerably lower scatters (36% lower for Eu, and 
6 1 % lower for Th). The [Th/Eu] abundance ratios thus obtained 
were used to determine the age of the Galactic disk in Paper II. 
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Fig. 12. Abundance patterns for all contaminating elements. Outliers are indicated by arrows labelled with their respective HD 
numbers. Solid lines, when present, represent exponential fits to the well behaved data points. 
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Fig. 14. Examples of spectral syntheses of the Eu ii line at 4129.72 A. FEROS and CES spectra are presented for two stars with 
extreme Eu abundances (HD 63 077 and HD 160691). Points are the observed spectra. Thick lines are the best fitting synthetic 
spectra, calculated with the shown Eu abundances. Thin lines represent variations in the Eu abundance A log e(Eu) = ±0.05 dex. 
The shown abundances are the ones that best fit the presented spectra, and not the average values that are presented in Table 1 141 
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Fig. 19. Examples of spectral syntheses of the Thn line at 4019.13 A. 3.60 m and CAT spectra are presented for two stars with 
extreme Th abundances (HD 20766 and HD 128 620). Points are the observed spectra. Thick lines are the best fitting synthetic 
spectra, calculated with the shown Th abundances. Thin lines represent variations in the Th abundance A log e(Th) = ±0.20 dex. 
The shown abundances are the ones that best fit the presented spectra, and not the average values that are presented in Tablell4l 



